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1.  Objectives 


The  main  objective  of  this  study  is  to  examine  coherent  quantum  effects,  such  as 
Rabi  oscillations  and  quantum  entanglement  in  optical  spectra  of  wide-band-gap 
materials,  and  to  determine  the  feasibility  of  fast  optical  control  of  quantum  states 
in  gallium  nitride  (GaN)  and  zinc  oxide  (ZnO)  heterostructures.  Because  of  much 
stronger  exciton-photon  interaction  in  these  materials  as  compared  to  gallium 
arsenide  (GaAs),  it  is  possible  to  realize  the  strong  coupling  regime  at  room 
temperature.  We  perform  an  experimental  and  theoretical  study  on  the  dephasing 
processes  in  available  nitride  and  zinc  oxide  materials.  This  will  allow  us  to  take 
advantage  of  their  unique  properties  for  observing  coherent  quantum  effects  at 
temperatures  higher  than  the  cryogenic,  at  liquid  nitrogen  temperature  and  higher. 
This  is  a  necessary  work  toward  coherent  optical  control  of  quantum  states  at 
higher  temperatures,  with  ultimately  room-temperature  coherent  control. 

1.1  Approach  to  Achieving  the  Objectives 

We  used  time-resolved  photoluminescence  (TRPL)  to  study  the  evolution  of 
carrier  density  and  interband  polarization  in  GaN  and  ZnO  materials.  Our 
modeling  of  ultrafast  optical  response  from  semiconductor  structures  is  based  on 
the  density  matrix  formalism.  The  detailed  numerical  solution  requires  a 
generalized  Monte  Carlo  simulation.  An  important  step  toward  coherent  control 
of  electronic  excitations  in  GaN  heterostructures  will  be  an  observation  of  Rabi 
oscillations  in  the  differential  transmission  signal.  In  the  wide-band-gap  materials 
like  GaN  and  ZnO,  the  exciton  transition  oscillator  strength  is  at  least  an  order  of 
magnitude  larger  than  in  GaAs;  therefore,  the  strong  coupling  regime  can  be 
achieved  at  significantly  higher  temperatures. 

1.2  Progress  Made  toward  Achieving  the  Director's  Research 
Initiative  (DRI)  Research  Objectives 

We  performed  theoretical  and  experimental  studies  of  the  sub-picosecond  kinetics 
of  photoexcited  carriers  in  GaN  and  ZnO.  In  the  theoretical  model,  interaction 
with  an  external  ultrafast  laser  pulse  is  treated  coherently,  and  to  account  for  the 
scattering  mechanisms  and  dephasing  processes,  a  generalized  Monte  Carlo 
simulation  is  used.  The  scattering  mechanisms  included  are  carrier  interactions 
with  polar  optical  phonons  and  acoustic  phonons,  and  carrier-carrier  Coulomb 
interactions.  We  studied  the  effects  of  various  scattering  mechanisms  on  the 
carrier  densities  and  also  presented  the  temperature  and  pulse-power-dependent 
luminescence  spectra.  The  results  are  presented  over  a  range  of  temperature. 


1 


electric  field,  and  excitation  energy  of  the  laser  pulse.  Here  we  report  the 
experimental  time-resolved  photoluminescence  studies  of  GaN  and  ZnO  samples. 
We  also  explain  the  kinetics  of  the  photoexcited  carriers  by  including  only  carrier- 
carrier  and  carrier-phonon  interactions  and  a  relatively  simple  2-band  electronic 
structure  model. 

The  analysis  and  modeling  of  polarization  evolution  and  coherent  oscillations  in 
GaN  quantum  wells  will  require  an  inclusion  of  carrier  confinement  effects  in  the 
quantum  kinetic  model.  Toward  this  goal,  we  have  included  polarization  fields 
and  strain  in  the  “k.p  model”  of  electron  and  hole  states  in  GaN  quantum  wells  as 
part  of  our  theory  of  radiative  recombination  rates  in  nitride  heterostructures.  The 
quantum  confinement  of  electron  states  is  required  for  effective  optical  control. 
Toward  that  goal,  we  modeled  the  evolution  of  the  spin  of  quantum-confined 
electrons,  providing  a  clear  geometrical  interpretation  of  the  qubit  dynamics. 
Using  the  analytic  form  of  the  evolution  operator,  we  proposed  a  set  of  single¬ 
qubit  rotations  that  is  solely  based  on  the  geometrical  phase.  We  demonstrated  the 
adiabatic  manipulation  of  a  qubit  using  only  the  geometric  phase.  This  has  some 
advantages,  since  it  reduces  the  requirements  of  perfect  tuning  of  the  controlled 
parameters  and  is  significantly  more  robust  against  noise. 

2.  Temperature  Effects  in  the  Kinetics  of  Photoexcited  Carriers 
in  Wide-Band-Gap  Semiconductors 

The  kinetics  of  photoexcited  electrons  and  holes  in  GaN  and  ZnO  were  studied  by 
sub-picosecond  time-resolved  photoluminescence.  The  rise  time  of  the  carrier 
recombination  was  analyzed  as  a  function  of  sample  temperature,  from  14  K  to 
room  temperature,  and  the  excitation  power.  In  one  set  of  measurements,  the 
excitation  wavelength  was  chosen  such  that  the  excess  photon  energy  was  750  meV 
above  the  band  gap,  well  above  the  polar  optical  phonon  energy.  In  another  set  of 
measurements,  the  excess  photon  energy  was  50  meV,  below  the  polar  optical 
phonon  energy.  The  goal  was  to  gain  a  better  understanding  of  the  scattering  and 
dephasing  processes  that  accompany  the  initial  photoexcitation  of  carriers.  Results 
showed  little  dependence  in  the  rise  time  on  the  initial  carrier  density  but  strong 
dependence  on  the  sample  temperature,  with  the  rise  time  increasing  by  more  than  1 
order  of  magnitude  as  temperature  decreases  from  room  temperature  down  to  14  K. 
To  elucidate  the  findings,  simulations  were  performed  based  on  a  2-band  model  and 
using  a  Monte  Carlo  technique.  In  the  model,  interaction  with  an  external  ultrafast 
laser  pulse  is  treated  coherently,  and  to  account  for  the  scattering  mechanisms  and 
dephasing  processes,  a  generalized  Monte  Carlo  simulation  is  used.  The  scattering 
mechanisms  included  are  carrier  interactions  with  polar  optical  phonons  and 
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acoustic  phonons,  and  carrier-carrier  Coulomb  interactions.  We  find  reasonable 
agreement  between  the  simulation  and  experiment. 

2.1  Theoretical  Model 

In  this  work  we  consider  a  bulk  wide-band-gap  semiconductor  material  described 
by  an  isotropic  parabolic  heavy  hole  band  and  an  isotropic  nonparabolic 
conduction  band.^  Considering  the  nonparabolicity  factor  for  the  conduction  band 
increases  the  accuracy  of  the  band  structure.  For  the  valence  band,  we  just 
consider  the  heavy  hole  band.  Within  the  energy  regions  considered  in  this  work, 
holes  will  be  confined  to  the  heavy  hole  band  because  of  its  higher  effective  mass, 
which  results  in  a  higher  density  of  states  compared  to  that  of  the  light  hole  band. 
Since  the  excitation  energies  that  we  consider  in  this  work  are  lower  than  the 
secondary  minima,  as  shown  in  Fig.  1,  the  particles  will  not  be  accelerated  or 
excited  to  energy  states  outside  the  F  valley. 
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Fig.  1  Constant  energy  surfaces  of  (a)  GaN  conduction  band  at  an  energy  of  750  meV  and 
(b)  GaN  valence  band  at  an  energy  of  100  meV  of  GaN  at  300  K 

Figure  la  shows  the  energy  surfaces  for  the  conduction  band  in  GaN  at  an  energy 
of  750  meV  above  the  conduction  band  minima,  and  Fig.  lb  shows  the  energy 
surfaces  for  the  valence  band  in  GaN  at  100  meV  below  the  valence  band 
maxima,  calculated  using  an  empirical  pseudopotential  model.  From  these  figures 
we  can  conclude  that  this  simple  2-band  model  is  a  good  approximation  of  the 
realistic  electronic  structure  for  the  range  of  energies  considered  in  this  work. 
However,  the  method  can  be  extended  to  multiple  valleys  or  full  band  structures. 
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The  Hamiltonian  describing  the  system  contains  2  parts:  Ho,  which  includes  the 
parts  that  can  be  treated  exactly  using  a  single  particle  model,  and  Hi,  which 
contains  various  interactions  and  dephasing  processes  that  will  be  treated  with 
some  level  of  approximations. 

In  this  work  we  consider  a  number  of  dephasing  processes  that  are  due  to  the 
interaction  of  carriers  and  interband  polarization  with  polar  optical  and  acoustic 
phonons,  and  the  carrier-carrier  interactions  via  Coulomb  potential.  We  describe 
the  dynamics  of  the  system  using  the  density  matrix  method.^’^  The  single  particle 
Hamiltonian  that  describes  the  interaction  of  free  carriers  with  the  incoming  light 
field  and  free  phonons  is  given  by 

k  k  k 

k 

where  fg  (1  +  ace)  =  k^/(2me)  and  Ch  =  h  /{Inih  )  are  the  non-parabolic  and 

parabolic  energy  states  for  electrons  and  holes,  and  me  and  m*  are  electron  and 
hole  effective  masses,  respectively.  Operators  cC,  dk^,(ck,  dk)  are  creation 
(annihilation)  operators  of  electrons,  holes  with  Fermi  anticommutation  relations. 
Operators  bk  and  bk^  are  phonon  operators  with  Bose  commutation  relations,  coq  is 
the  phonon  frequency,  Mk  is  the  interband  dipole  matrix  element,  and  £"0(0  is  the 
amplitude  of  the  Gaussian-shaped  laser  pulse  with  a  frequency  centered  at  (Ol  and 
defined  as  £0  (t)  =  El  exp(-f/T  ^),  where  El  is  the  peak  electric  field  and  z  can  be 
related  to  the  full  width  at  half  maximum  of  the  measured  field  intensity  by  r  =  ri 
/(2^2ln2).  The  external  light  field  is  treated  semi-classically  within  the  rotating 
wave  approximation."^  The  optical  dipole  matrix  element  in  the  parabolic  bands 
approximation  is  given  by  Mk  =  dcv(I+  h^k^/lmrEg)'^ ,  which  can  be  approximated 
by  a  constant,  Mk  ~  dev,  for  large  band  gap  values  Eg.  The  reduced  mass  is  defined 
by  mr  =  (me  ^  -i-  mh^yK  The  interband  electric  dipole  dev  is  related  to  interband 
momentum  matrix  element  dev/e  =  hpcv/Egmo  and  in  the  effective  mass 
approximation  for  2-band  model  dev  =  [(-l+mo/me)/2Egmo)]^^,  from  which  we 
obtain  in  GaN  case  dev/e  2.12  A. 


To  describe  the  kinetics  of  the  system,  we  need  to  obtain  the  distribution  function 
of  the  carriers  and  polarization  as  well  as  the  rate  of  change  of  these  functions. 


d  - 

Using  the  Heisenberg  equations  of  motion,  ifi  —  0{t)  = 

dt 


0,H 


for  the  electron  and 


hole  operators  and  incorporating  appropriate  commutation  relations,  we  obtain  the 
equation  of  motion  for  particles  and  polarization  that  represent  the  generation  rate 
for  a  noninteracting  system: 
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(2) 

(3) 

(4) 


where  fk,  fk  and  pk  are  the  distribution  functions  of  electrons,  holes,  and 
polarization,  respectively: 

fk  it)  =  (c;  {t)c,  (O),  A"  it)  =  {d:  (t)d,  (t)),  p,  (t)  =  {d_,  (t)c,  (0) ,  with  brackets 
denoting  the  expectation  value. 

Various  interactions  among  the  particles  are  introduced  in  the  model  through  a 
perturbing  Hamiltonian.  For  each  interaction  mechanism,  we  consider  a 
Hamiltonian  Hf  where  i  represents  the  specific  interaction  mechanism.  In  this 
work  we  consider  polar  optical  and  acoustic  phonons,  and  carrier-carrier 
interactions.  A  number  of  suitable  approximations  are  also  included  to  overcome 
the  difficulties  inherent  in  the  many-body  nature  of  these  interactions.  The 
carrier-carrier  interaction  Hamiltonian  is  given  by 


-c:d:^d_ 

k,PS  ^ 


p+q  k-q  •> 


where  the  first  and  second  terms  are  the  repulsive  electron-electron  and  hole-hole 
interactions,  and  the  third  term  is  the  attractive  interaction  of  electrons  and  holes, 
and  Vq  is  the  screened  Coulomb  potential.  The  first-order  approximation  in  the 
Heisenberg  equations  of  motion  leads  to  a  renormalization  of  the  carrier  energies 
by  a  self-energy  factor  and  renormalization  of  the  light  field  by  an  internal  field.  ^ 
After  obtaining  the  first-order  correction,  using  the  Heisenberg  equation  of 
motion  again,  one  can  derive  the  equations  of  the  second-order  contribution  of 
carrier-carrier  interaction.  The  equation  for  expectation  values  of  4  operators 
involves  the  strings  of  6  operators.  We  use  factorizations  to  truncate  these  chains 
of  equations.^ 

The  Markov  approximation^  "^  is  used  to  avoid  time  integral  equations  and  to 
neglect  the  retardation  effects.  This  approximation  also  implies  that  the  electric 
field  and  the  distribution  functions  vary  slowly  compared  to  the  oscillation  of  the 
exponential  terms  of  the  laser  pulse;  therefore,  it  eliminates  the  effect  of  initial 
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correlations  and  dynamic  memory.  Furthermore,  the  validity  of  this  approach  is 
limited  to  the  range  of  carrier  densities  considered  in  this  work  and  may  need  to 
be  revised  for  higher  densities.  In  this  work  we  consider  the  second-order 
contribution  of  the  carrier-carrier  interaction.  The  result  of  the  second-order 
carrier-carrier  interaction  excluding  the  Hartree-Fock  terms  are  given  in  Eqs.  6 
and  7,  where  a  =  e,h: 


dfk 


a 


cc(2) 


dt 


a(cc) 


k-q.kfk 


(6) 


dPk 


cc{2) 


dt 


-Zk'*  ■9<kPi^  —  W^k,k-qPi^_g^, 


(7) 


where  the  carrier-carrier  transition  matrices  are  given  by 


af  kr 

-  fk-q)  [/fc“'(l  -  fk'+q)  -  Pl'+qPk]  +  C.C., 
and  the  polarization  transition  matrices  are  given  by 
^k-q,k  =  X  I"  Z  Z  ~  ^  [~P*k'+c,Pk' 

a,ar  k' 

+ fs'  (i  -  4-;,)  (1  -  ff-,) + -  //?')]  ■ 

The  function  T)  (c)  is  defined  as  follows; 

1  T 

D(oo)  =  : - 1-  5(g))  , 

m  G) 


(8) 


(9) 


(10) 


where  T  denotes  principal  value. 


Because  of  the  presence  of  free  carriers,  the  Coulomb  potential  is  screened.  We 
use  a  static  screening  potential 


c"  1 


(11) 


where  co  and  Cs  are  the  vacuum  permittivity  and  static  dielectric  constant, 
respectively.  The  screening  wave  vector  is  calculated  according  to  the  following 
equation: 
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where  e(k)  is  the  carrier  energy  state,  and  V  is  the  volume  of  the  material. 


Besides  the  interaction  of  carriers  with  each  other,  in  this  work  we  take  into 
account  the  interaction  with  polar  optical  and  acoustic  phonons  as  well.  The 
structure  of  the  Hamiltonian  describing  these  interactions  is  independent  of  the 
type  of  the  phonon  branch.  The  only  difference  is  the  interaction  matrix  element 
that  describes  the  type  of  lattice  deformation,  in  the  case  of  acoustic  phonons,  and 
the  long-range  electric  field  when  polar  optical  modes  are  considered.  The 
Hamiltonian  that  describes  the  carrier-phonon  interaction  is 


+r;’ctb;c,,,  +r;d;„b,d,  +r‘;dtb;d,^],  (13) 


where  q  is  the  interaction  matrix  element.  Because  of  the  opposite  charge  of 
electrons  and  holes,  the  matrix  elements  in  the  case  of  a  polar  interaction  are 
f  q  =  q  =  jq.  The  proccdurc  of  obtaining  the  second-order  contribution  of  the 
carrier-phonon  interaction  to  the  equations  of  motion  is  similar  to  that  of  carrier- 
carrier  interaction.  We  derive  Heisenberg  equations  of  motion  using  the  Markov 
approximation.  The  second-order  contribution  to  the  equations  of  motion, 
excluding  the  internal  field  terms,  is  the  same  as  Eqs.  6  and  7;  however,  instead  of 
Yi/e;h(cc)  ^/p(cc) ^  Contribution  of  phonons  and  must  be 

substituted.  The  carrier-phonon  scattering  rates  are  given  by 


,,,a,(cp)  _  ^ 
^\-q,k  -  ft 


I.±-\Yq\^K^k 


—  8 


a 

k-q 


±  hOJq)  {Nq+^±  (l  -  fk-q). 


(14) 


and  the  polarization  transition  matrices  are  given  by 


<•?  =  i  S  ±  [(«, + i 

a=e,h  ±- 


(15) 


where  Nq  is  the  number  of  phonons  with  wave  number  ^  at  a  specific  temperature, 
which  is  calculated  using  the  Bose-Einstein  distribution  and  coq  is  the  phonon 
frequency. 

To  obtain  the  distribution  functions,  the  equations  are  separated  into  2  parts,  a 
coherent  part  and  a  dephasing  part.  The  coherent  contribution  includes  the 
generation  equations  that  are  calculated  according  to  Eqs.  3  and  4  for  particles  and 
carriers,  respectively.  The  value  of  the  distribution  functions  in  the  coherent  case 
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can  be  obtained  by  directly  integrating  the  generation  equations.  A  generalized 
Monte  Carlo  method,^  which  is  very  similar  to  the  traditional  Monte  Carlo 
method,  is  used  to  account  for  the  dephasing  phenomena.  These  distribution 
functions  are  then  used  to  calculate  the  luminescence  spectrum.  The  rate  of  the 
scattering  mechanisms  employed  in  the  Monte  Carlo  procedure  is  evaluated  using 
the  Fermi  golden  rule.  One  of  the  most  important  scattering  mechanisms  is  the 
carrier-carrier  interaction,  especially  in  GaN  with  carrier  densities  as  high  as 
10^^  cm“^  as  considered  in  this  work.  At  this  range  of  carrier  densities,  the 
screening  effect  is  essential  both  for  the  carrier-carrier  scattering  and  carrier- 
phonon  scattering.  In  fact,  each  carrier,  while  interacting  with  other  carriers 
through  the  screened  Coulomb  potential,  is  also  influenced  by  the  screening 
potential  in  its  interactions  with  phonons.  Therefore,  the  carrier-phonon 
interaction  is  weakened  considerably  by  the  screening  potential.  The  screened 
carrier-carrier  scattering  is  calculated  based  on  the  model,^’^  which  takes  into 
account  the  static  screening  wave  vector  as  shown  in  Eq.  12.  The  matrix  element 
of  the  interaction  of  the  electrons  with  the  polar  optical  phonons  via  the  Frohlich 
interaction  is  as  follows^: 


r. 


e  ficOjjj 

r  1 

_  0 

1 

^0 

'(e+K-f 

-r{k,k  +  q)), 


(16) 


where  Coo  is  the  high-frequency  dielectric  constant,  q  is  the  wave  vector  of 
phonons,  cow  is  the  polar  optical  phonon  energy,  and  /(k,k’)  is  the  overlap  integral 
for  the  conduction  band.^  For  the  holes,  we  consider  the  ordinary  polar  optical 
scattering  rate  for  parabolic  band  structures.  The  calculated  polar  optical 
scattering  rate  at  room  temperature,  including  the  effect  of  screening,  is  shown  in 
Fig.  2. 
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Fig.  2  Calculated  polar  optical  scattering  rate  for  the  nonparabolic  conduction  band  in 
GaN  including  the  screening  effect.  The  solid  lines  show  the  phonon  emission  rate  and  the 
dashed  lines  show  the  phonon  absorption  rate.  The  electric  field  of  the  light  is  0.6  x  10^  Vm“^ 
and  the  excess  photon  energy  is  750  meV. 

The  excess  photon  energy — the  amount  of  photon  energy  after  subtracting  the 
band  gap  energy — is  750  meV,  and  the  electric  field  of  the  external  light  is  0.6  x 
10^  Vm“^  The  solid  lines  show  the  phonon  emission  rates,  and  the  dashed  lines 
show  the  rate  for  the  phonon  absorption  processes.  The  black  lines  (topmost  solid 
and  dashed  lines)  show  the  calculated  scattering  rate  at  -100  fs.  At  this  time,  the 
generation  process  has  just  started,  and  the  density  of  electrons  in  the  system  is 
8.2  X  10^^  cm"^.  The  blue  lines  (center  solid  and  dashed  lines)  show  the  scattering 
rates  at  0  fs  when  the  maximum  of  the  laser  pulse  occurs  and  the  density  of 
electrons  is  6.1  x  10^^  cm“^.  Because  of  the  increasing  carrier  density  resulting 
from  a  higher  generation  rate,  the  screening  wave  vector  changes  dynamically  as 
the  simulated  system  evolves.  Therefore,  the  polar  optical  scattering  will  be 
screened,  and  eventually  the  scattering  rate  will  decrease.  This  effect  is  more 
pronounced,  as  shown  by  the  red  lines  (bottom  solid  and  dashed  lines),  when  the 
time  is  100  fs  and  more  carriers  are  generated.  In  this  case,  the  density  of 
electrons  is  1.2  x  10^^  cm“^,  and  the  electron-phonon  scattering  rate  is  at  its 
lowest  value.  Without  the  screening  effect,  the  polar  optical  scattering  rate  will 
not  change  during  the  simulation  run. 

The  deformation  potential  acoustic  scattering  rate  for  isotropic  parabolic  and 
nonparabolic  bands  are  calculated  based  on  the  method  presented  in  Jacoboni  and 
Lugli.^  For  the  parabolic  bands,  only  longitudinal  modes  contribute  to  the 
scattering.  The  scattering  rate  is  calculated  considering  the  acoustic  deformation 
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potentials  {De  and  Dh),  longitudinal  sound  velocity  (Vt),  density  of  the  material 
(p),  and  the  phonon  number,  nq,  which  is  calculated  based  on  Bose-Einstein 
distribution.  For  the  nonparabolic  conduction  band,  the  average  value  of  the 
longitudinal  and  transverse  sound  velocity  {Vt  )  is  used.  The  transferred  energy  is 
calculated  based  on  energy  and  momentum  conservation.  The  material  parameters 
that  have  been  used  in  this  work  are  listed  in  Shishehchi  et  al.^’^ 

2.2  Numerical  Procedure 

This  section  describes  the  numerical  approach^  used  to  evaluate  the  distribution 
functions  based  on  the  model  presented  in  Section  2.1.  Based  on  this  model,  the 
ensemble  of  the  carriers  is  described  by  means  of  a  number  representation.  This 
means  that  the  number  of  super  particles  representing  the  carriers  is  defined  at 
each  cell  in  k  space;  therefore,  the  occupation  of  each  cell  is  decreased  or 
increased  by  one,  depending  on  the  mechanism  that  takes  place.^  Furthermore,  the 
maximum  occupancy  of  each  cell  is  limited  by  the  density  of  states.  The 
distribution  function  is  then  defined  by  dividing  the  occupation  number  of  each 
cell  by  the  density  of  state.  The  distribution  functions  of  electrons  and  holes  are 
real  variables,  whereas  that  of  the  polarization  is  complex.  Consequently,  the 
magnitude  and  the  phase  of  polarization  are  assigned  to  a  2-vector  matrix.  The 
simulation  process  initiates  with  the  vacuum  state  of  the  distribution  of  carriers 
and  polarization;  in  other  words,  the  distribution  functions  and  polarization  are 
equal  to  zero.  Subsequently,  the  ensemble  of  the  carriers  and  polarization  evolves 
as  time  elapses  as  a  result  of  the  laser  pulse  taking  the  system  out  of  equilibrium. 

The  simulation  period  is  divided  into  a  number  of  time  steps.  At  each  time  step, 
first  the  coherent  part  is  calculated.  The  generation  rate  for  the  carriers  is 
evaluated  by  a  direct  anal3^ical  integration,  and  the  time  evolution  of  the 
polarization  is  obtained  by  integrating  the  coupled  differential  Eqs.  2  and  3 
employing  a  Range-Kutta  algorithm.  The  obtained  distributions  until  this  point 
are  then  transformed  into  a  many-body  particle  distribution,  which  can  be 
efficiently  used  to  implement  the  ensemble  Monte-Carlo  calculation.  Specifically, 
we  use  N  super-particles  for  each  kinetic  variable.  These  particles  are  generated 
randomly  according  to  the  distribution  function  calculated  during  the  coherent 
part  of  the  time  step.  Since  the  polarization  distribution  function  is  complex,  the 
phase  is  discretized  according  to  the  complex  values  of  pk.  Next,  the  scattering 
rates  are  calculated.  Because  of  the  dependence  of  the  screening  length  on  the 
distribution  functions,  which  are  updated  at  each  time  step,  scattering  rates  should 
be  recalculated  once  the  new  distributions  are  available.  Then,  a  traditional  Monte 
Carlo  simulation  is  performed  using  the  newly  calculated  scattering  rates.  After 
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each  scattering  event,  the  number  of  super  particles  in  the  initial  cell  is  decreased 
by  one,  and  the  number  in  the  final  cell  is  increased  by  one. 

In  this  work,  we  also  include  the  effect  of  the  Pauli  exclusion  principle.  Although 
the  Monte  Carlo  method  treats  the  electrons  and  holes  as  semi-classical  particles, 
they  should  obey  the  Pauli  exclusion  principle,  especially  at  high  carrier  densities 
where  the  distribution  of  the  particles  is  highly  influenced  by  their  fermionic 
properties. After  the  scattering  process,  once  the  final  state  is  defined,  a  random 
number  between  0  and  1  is  generated  to  evaluate  whether  the  transition  is 
accepted  or  rejected. In  this  method,  the  distribution  function,/; ,  at  each  cell,  i, 
in  the  k  space  is  compared  to  the  random  number,  r.  If  r  >  fi,  the  transition  is 
accepted  and  the  occupancy  of  that  cell  is  increased  by  one.  Otherwise,  the 
transition  is  denied  and  the  particle  goes  through  a  self-scattering  event. 
Eventually,  the  distribution  function  is  updated  according  to  the  exclusion 
principle  after  scattering  by  each  mechanism. 

2.3  Experiment 

A  schematic  representation  of  the  experimental  setup  is  shown  in  Fig.  3.  A 
250-kHz  regeneratively  amplified  Ti:Sapphire  laser  produces  150-fs  pulses  at 
800-nm  wavelength.  Part  of  this  beam  is  split  off,  acting  as  a  gating  pulse,  and  the 
other  part  pumps  an  optical  parametric  amplifier  tunable  from  760  to  465  nm. 
Output  from  the  amplifier  is  again  doubled  in  a  nonlinear  optical  crystal  (2Xtal)  to 
produce  pulses  tunable  from  380  to  232.5  nm.  The  pump  pulse,  compressed  to  a 
width  of  approximately  80  fs,  is  focused  onto  the  sample  where  it  photo-excites 
electron-hole  pairs.  The  luminescence  from  the  sample  is  collected  using 
nondispersive  optics  and  focused  on  another  nonlinear  crystal  (NLXtal),  where  it 
is  mixed  with  the  gating  pulse  to  yield  a  down-converted  signal  that  is  spatially 
and  spectrally  separated  from  the  other  beam.  The  down-converted  signal  is 
passed  through  a  spectrometer  and  detected  with  a  photomultiplier  tube  in  photon 
counting  mode.  The  sample  used  in  this  experiment  was  a  2-pm-thick 
heteroepitaxial  GaN  film  deposited  by  metalorganic  chemical  vapour  deposition 
(MOCVD)  with  a  fully  graded  aluminum  gallium  nitride  (AlGaN)  layer  on  bulk 
aluminum  nitride  (AIN).  The  arrival  time  of  the  near-IR  gate  pulse  relative  to  the 
UV  photoluminescence  maps  the  TRPL  with  350-fs  resolution  as  determined  by 
the  group  velocity  dispersion  between  the  2  beams  in  the  nonlinear  crystal.  The 
energy  range  of  the  carrier  distribution  monitored  is  a  function  of  the  acceptance 
angle  of  the  phase  matching  condition  in  the  nonlinear  crystal. 
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Fig.  3  A  schematic  representation  for  time -resolved  photoluminescence  through  optical 
gating  by  frequency  down-conversion.  Excitation  laser  provides  pulses  tunable  from  3.26  to 
5.33  eV,  suitable  for  studying  GaN  with  generated  electron-hole  pairs  up  to  10^®  cm“^. 

2.4  Results  and  Discussion  for  the  Gallium  Nitride  (GaN)  Case 

In  this  section  we  present  and  compare  the  theoretical  and  the  experimental  results 
for  the  case  of  GaN.  In  particular,  we  will  look  at  the  time  evolution  of  the 
electron  and  hole  densities  as  a  function  of  the  electric  field  strength  and  intensity 
of  the  external  pulse.  The  material  parameters  are  listed  in  Table  1.  We  present 
the  effect  of  various  scattering  mechanisms  on  these  results,  with  particular 
emphasis  on  the  analysis  of  the  phenomena  that  lead  to  the  disruption  of  the 
coherence  in  the  system.  In  this  context,  to  highlight  the  role  of  polar  optical 
phonon  scattering,  we  consider  2  excitation  energies:  one  with  an  excess  photon 
energy  of  750  meV,  which  is  higher  than  the  polar  optical  phonon  energy  (91.2 
meV),  and  the  other  one  with  excess  photon  energy  of  50  meV,  which  is  lower 
than  polar  optical  phonon  energy.  For  the  case  of  excitation  with  the  lower 
energy,  longitudinal  optical  (LO)  phonon  emission  processes  are  significantly 
suppressed  and  do  not  drive  the  energy  relaxation  in  the  system.  Finally,  we 
consider  the  calculated  luminescence  for  both  values  of  the  excess  photon 
energies.  Moreover,  we  study  the  effect  of  temperature  on  the  luminescence 
spectra  and  compare  it  with  the  measured  values. 
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Table  1  Parameters  used  in  the  simulation  for  GaN 


Rarameter 

Valiie 

Unit 

0.22inei 

kg 

kg 

too 

K  m  ’ 

F  m  ^ 

hfiHQ 

91  2 

nieV 

3.39 

eV 

if 

2.12 

A 

a 

0.37 

eV  ’ 

D, 

H.3 

eV 

22 

eV 

P 

6  15 

gcni  ^ 

Vl 

7.96 

lO^cnis  ’ 

Vj 

4  13 

lO^cnis  ’ 

T 

}^5 

fs 

Figure  4a  shows  the  calculated  time  evolution  of  the  electron  number  density 
function  at  room  temperature  when  the  excess  photon  energy  is  750  meV  and  the 
electric  field  of  the  laser  pulse  is  0.6  x  10^  V  m“^  The  number  density  function 
n{EI)  is  defined  as  the  product  of  the  density  of  states  and  the  carrier  distribution: 
n{E)  =  g{E)  f(E).  Figure  4b  provides  the  same  information  for  holes.  The  energy 
in  (a)  is  measured  from  the  bottom  of  the  conduction  band,  and  the  energy  in  (b) 
is  measured  from  the  top  of  the  valence  band. 


Fig.  4  Calculated  time  evolution  of  the  electron  (a)  and  hole  (b)  number  density  functions 
at  room  temperature 
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The  time  scale  has  been  chosen  so  that  the  peak  of  the  Gaussian  laser  pulse  occurs 
at  0  fs,  and  the  calculated  electric  field  interaction  extends  from  -200  to  200  fs. 
The  particle  distribution  results  are  shown  at  the  instants  corresponding  to  -50  fs, 
0  fs,  50  fs,  100  fs,  150  fs,  and  500  fs.  The  dashed  black  line  shows  the  normalized 
electron  density  at  -50  fs.  At  this  time  step,  carriers  are  being  generated  and,  as 
the  time  elapses,  the  number  of  the  carriers’  density  increases.  Since  the  excess 
photon  energy  of  750  meV  is  divided  between  the  conduction  and  valence  bands 
(based  on  the  respective  density  of  states),  we  observe  that  the  maximum  of  the 
electron  distribution  is  below  the  excess  photon  energy.  For  holes,  the  peak  of  the 
distribution  occurs  just  below  200  meV.  As  the  time  elapses,  the  scattering 
processes  relax  the  carriers’  energy  and  change  the  shape  of  the  distribution.  The 
dot-dashed  green  line  shows  the  results  at  0  fs — the  time  at  which  the  laser  pulse 
intensity  is  maximum.  At  this  time,  another  peak  emerges  at  an  energy  of  about 
91  meV,  which  corresponds  to  the  optical  phonon  energy.  The  peak  is  due  to  the 
polar  optical  phonon  scattering  emission  mechanism  that  cannot  occur  below  this 
energy.  At  later  times,  100  and  150  fs,  the  shoulder  at  =^91  meV  becomes  more 
and  more  pronounced.  Eventually,  at  500  fs  when  the  laser  pulse  is  extinguished, 
the  distribution  is  completely  relaxed  and  carriers  have  accumulated  in  a  region 
close  to  the  band  edge  energy. 

From  the  time  evolution  of  the  carrier  number  density  functions,  it  is  possible  to 
gamer  additional  information  on  the  coherent  processes  in  the  system.  The 
coherent  effects  can  be  analyzed  either  as  a  function  of  the  electric  field  at  a  given 
time  or  for  a  given  electric  field  intensity  as  a  function  of  time.  Figures  5a  and  5b 
show  the  calculated  normalized  total  electron  density  as  a  function  of  the  electric 
field  of  the  external  light  for  the  case  that  the  excess  photon  energy  is  50  and  750 
meV,  respectively.  In  this  case,  the  lattice  temperature  is  300  K,  and  the  density 
values  are  obtained  after  the  system  has  reached  steady  state,  long  after  the  laser 
pulse  is  finished  (1.5  ps).  Furthermore,  the  integrated  number  of  the  electrons  in 
the  whole  k  space  has  been  calculated,  then  the  data  have  been  normalized  to  the 
maximum  value  obtained  at  each  field  value. 
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Fig.  5  (a)  Calculated  effect  of  scattering  mechanisms  on  the  integrated  density  of 

electrons  in  GaN  vs.  electric  field  of  the  external  light  at  300  K  when  the  excess  photon 
energy  is  50  meV.  (b)  Calculated  effect  of  scattering  mechanisms  on  the  integrated  density  of 
electrons  in  GaN  vs.  electric  field  of  the  pulse  at  300  K  when  the  excess  photon  energy  is 
750  meV. 

The  various  curves  in  Figs.  5a  and  5b  illustrate  the  effect  of  various  scattering 
mechanisms  on  the  density  of  the  carriers  and  degree  of  coherence  in  the  system. 
The  dashed  black  line  shows  the  coherent  case  where  no  dephasing  process  is 
included.  In  this  case,  the  Monte  Carlo  part  of  the  numerical  model  is  excluded, 
and  the  results  are  obtained  by  just  taking  into  account  the  first-order  solution  of 
the  equations  of  motion.  Since  no  dephasing  process  is  present,  the  distribution 
functions  of  electrons  and  holes  are  highly  correlated.  Furthermore,  Rabi 
oscillations"^  can  be  seen  in  the  result.  The  oscillation  is  more  pronounced  when 
the  excitation  energy  is  higher  (750  meV  in  Fig.  5b).  Figure  6  presents  the  Rabi 
oscillation  in  the  time  domain  for  the  2  excess  energy  values. 
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Fig.  6  Calculated  normalized  electron  density  vs.  time  at  room  temperature  in  GaN.  The 
electric  field  of  the  light  pulse  is  10*  Vm"^  for  2  values  of  excess  photon  energy.  These  curves 
have  been  computed  only  considering  the  coherent  effects  as  described  in  Eqs.  2,  3,  and  4, 
and  neglecting  the  dephasing  due  to  scattering  mechanisms.  The  values  are  normalized  to 
the  maximum  of  the  curve  corresponding  to  750  meV. 

The  normalized  electron  density  at  room  temperature  as  a  function  of  time  for  the 
previous  excitation  energies  when  the  electric  field  is  equal  to  10^  Vm"^  is  shown 
in  Fig.  6.  The  excess  photon  energies  for  the  solid  and  dashed  lines  are  750  and  50 
meV,  respectively.  It  can  be  seen  that  the  Rabi  oscillation  becomes  stronger  as  the 
excitation  energy  increases.  Besides  the  case  in  which  no  dephasing  is  present. 
Fig.  5  shows  the  effects  of  various  scattering  processes  on  the  coherence  in  the 
system.  The  presence  of  interactions  with  phonons  and  other  carriers  changes  the 
picture  of  the  correlation  between  electrons  and  holes.  The  dot-dashed  green  lines 
in  Figs.  5a  and  5b  show  the  case  that  only  phonon-scattering  mechanisms  are 
included.  The  phonon-scattering  process  includes  both  the  polar  optical  and 
acoustic  scattering  mechanisms.  In  polar  optical  phonon  interaction,  the  carriers 
emit  or  absorb  an  integer  number  of  phonons,  and  the  distribution  relaxes  toward 
lower  energy  states.  Consequently,  if  the  excitation  energy  is  lower  than  the  polar 
optical  energy,  the  polar  optical  phonon  emission  process  is  less  probable. 
However,  acoustic  phonon  processes,  having  a  much  lower  energy  than  polar 
optical  phonons,  can  occur  at  any  excitation  energy.  The  dashed  blue  lines 
represent  the  case  when  only  the  carrier-carrier  interaction  is  included  in  the 
system.  The  role  of  this  process  is  to  basically  redistribute  the  energy  in  the 
system;  therefore,  as  a  result  of  this  relaxation  mechanism,  the  peaks  of  the  carrier 
distribution  function  broaden.  Eventually,  the  solid  red  line  represents  the  case 
that  all  the  above  mechanisms  are  included.  The  broadening  effect  due  to  the 
carrier-carrier  scattering  is  dominant  when  the  excess  photon  energy  is  50  meV 
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since  the  polar  optical  is  no  longer  the  dominant  energy  loss  mechanism.  To 
summarize,  according  to  Fig.  5a,  for  the  case  that  the  excess  photon  energy  is  less 
than  the  polar  optical  phonon  energy,  since  the  polar  optical  phonon  emission  is 
no  longer  dominant,  the  total  density  of  electrons,  where  all  the  interactions  are 
present,  follows  the  shape  of  the  carrier-carrier  dephasing  curve.  This  suggests 
that  in  this  case,  the  carrier-carrier  scattering  mechanism  is  the  dominant  process. 

However,  in  Fig.  5b  where  the  photon  excess  energy  is  higher  than  the  polar 
optical  phonon  energy,  polar  optical  phonon  emission  is  a  strong  process  in 
relaxing  the  carriers;  therefore,  the  final  density  of  electrons  follows  the  phonon 
dephasing  curve.  In  this  figure,  the  electron  density  values  of  the  carrier-carrier 
dephasing  case  at  each  electric  field  seem  to  be  larger  than  the  values  of  the  curve 
where  all  the  scattering  rates  are  included.  However,  care  must  be  taken  in 
comparison  since  these  are  the  normalized  values.  The  actual  values  of  the  carrier 
density  at  each  electric  field  in  the  case  where  all  the  scattering  mechanisms  are 
included  exceed  the  values  of  all  the  other  cases.  In  fact,  the  presence  of  each 
scattering  mechanism  affects  the  distribution  of  polarization,  which  has  a  major 
role  in  the  carrier  generation  rate. 

To  compare  the  results  obtained  from  the  numerical  model  with  the  experimental 
data,  we  have  computed  the  energy-resolved  time-dependent  luminescence.  To 
calculate  the  luminescence  at  various  energies  of  interest,  once  the  system  reaches 
the  steady  state  condition,  we  select  specific  ranges  of  photon  energies.  For 
example,  we  consider  narrow  energy  regions  close  to  the  band  gap  energy,  the 
polar  optical  phonon  energy,  twice  the  polar  optical  phonon  energy,  and  other 
multiples.  Subsequently,  we  multiply  the  electron  and  holes  at  each  energy  within 
the  narrow  band  and  sum  all  these  partial  products  over  the  selected  energy 
region.  After  that,  we  normalize  the  luminescence  data  to  the  respective 
maximum  values.  Figure  7  shows  the  calculated  energy-resolved  time-dependent 
photoluminescence  for  the  case  that  the  excess  photon  energy  is  750  meV.  The 
electric  field  of  the  external  pulse  is  0.6  x  10^  Vm"^  and  the  lattice  temperature  is 
300  K.  The  peak  of  the  laser  pulse  is  located  at  the  origin  of  the  time  axis,  and  the 
time  evolution  of  the  carrier  ensembles  has  been  collected  up  to  8  ps. 
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Fig.  7  Luminescence  spectra  in  GaN  at  room  temperature  (300  K)  calculated  for  different 
energy  regions 

The  inset  of  the  figure  shows  the  results  up  to  0.6  ps  where  the  details  and  the 
sequence  of  the  curves  are  more  distinguishable.  The  solid  black  line  shows  the 
results  when  the  photoluminescence  energy  is  collected  at  band  edge  energy, 
which  is  the  slowest  process.  The  dotted  red  line,  the  dot-dashed  green  line,  and 
the  dashed  blue  line  show  the  results  when  the  selected  energy  is  equal  to  1,  2, 
and  3  polar  optical  phonon  energies  (ILO,  2LO,  3LO),  respectively.  The  3LO 
case  is  the  fastest  process,  then  the  2LO,  and  afterwards  the  ILO  curves  reach  the 
steady  state. 

Figures  8a  and  8b  present  the  luminescence  results  in  GaN  when  the  temperature 
is  equal  to  50  and  150  K,  respectively.  The  electric  field  and  the  excess  photon 
energy  are  the  same  as  the  ones  in  Fig.  7. 

The  results  clearly  show  that  as  the  temperature  decreases,  the  rise  time  increases. 
More  specifically,  the  rise  time  for  the  band  edge  energy  region  changes 
significantly.  This  is  due  to  the  fact  that  at  lower  temperatures,  the  polar  optical 
scattering  rate  is  lower,  and  therefore  this  process  is  weaker  in  relaxing  the  carrier 
distributions. 

Although  we  have  endeavored  to  develop  our  simulation  models  to  provide  the 
best  description  of  the  physical  phenomena  that  emerge  from  the  experiment,  a 
quantitative  comparison  is  not  always  possible.  In  fact,  the  simplified  band 
structure  and  neglecting  the  internal  field  and  the  imaginary  part  of  the  self¬ 
energy  arising  from  the  second-order  correction  to  the  carrier-carrier  and  carrier- 
phonon  interactions^  may  lead  to  simulation  results  that  in  some  cases  can  only  be 
in  qualitative  agreement  with  the  experiment. 
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Fig.  8  (a)  Luminescence  spectra  in  GaN  at  150  K  collected  at  different  energy  regions. 

The  excess  photon  energy  is  750  meV.  (b)  Luminescence  spectra  in  GaN  at  50  K  collected  at 
different  energy  regions.  The  excess  photon  energy  is  750  meV.  The  electric  field  of  the 
external  light  pulse  is  0.6  x  10’  V  m“^  in  both  cases. 

Figures  9a  and  9b  show  the  results  of  luminescence  for  the  temperature  values 
ranging  from  150  to  300  K  obtained  by  experiment  and  simulation,  respectively, 
when  the  excess  photon  energy  is  750  meV.  It  clearly  shows  that  as  the  lattice 
temperature  increases,  the  luminescence  rise  time  decreases. 

Because  of  the  specific  experimental  conditions,  it  is  not  always  possible  to 
directly  correlate  the  beam  power  used  in  the  measurements  with  the  electric  field 
magnitude  employed  in  the  simulations.  In  fact,  for  the  model  we  have  employed, 
as  opposed  to  the  case  of  a  semiclassical  generation  model,  the  number  of 
generated  carriers  depends  on  the  interplay  between  the  coherent  and  dephasing 
phenomena.  Nevertheless,  the  numerical  model  predicts  the  correct  trend  of  the 
experimental  results.  It  is  important  to  appreciate  this  point,  since  the 
experimental  conditions  under  which  the  results  in  Fig.  9a  were  obtained  are  not 
the  same  as  the  ones  considered  in  computing  the  values  presented  in  Fig.  9b. 
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Fig.  9  Luminescence  spectra  of  GaN  at  different  temperatures  collected  at  band  edge 
energy  region.  The  excess  photon  energy  is  750  meV.  (a)  shows  the  experimental  results  and  (b) 
presents  the  simulated  results  for  an  electric  field  of  the  external  light  pulse  of  0.6  x  10^  V  m"^- 


In  order  to  show  that  this  is  still  a  valid  comparison,  even  if  the  power  of  the  beam 
in  the  experiment  and  the  electric  field  magnitude  in  the  numerical  model  cannot 
be  directly  correlated,  we  have  computed  the  luminescence  as  a  function  of  the 
electric  field  strength.  According  to  the  experimental  results  shown  in  Fig.  10a, 
the  luminescence  rise  time  is  independent  of  the  power  of  the  laser  pulse.  In  both 
Figs.  10a  and  10b,  the  luminescence  results  at  GaN  collected  at  the  band  edge 
energy  region  are  shown  over  a  range  of  laser  pulse  powers  in  which  the  rise 
times  are  approximately  the  same. 
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Fig.  10  Luminescence  spectra  of  GaN  at  different  laser  powers  collected  at  band  edge 
energy  region.  The  excess  photon  energy  is  750  meV  and  the  lattice  temperature  is  150  K.  (a) 
shows  the  experimental  results  and  (b)  presents  the  simulated  results.  In  (a),  E  represents 
the  baseline  power  used  in  the  experiment.  In  (b),  E  is  the  baseline  value  of  electric  field 
employed  in  the  simulation. 

In  Fig.  10,  the  lattice  temperature  is  150  K  and  the  excess  photon  energy  is 
750  meV.  Fig.  10b  shows  that  the  numerical  simulation  results  deliver  the  same 
trend  as  the  experiment.  On  the  basis  of  this  result,  we  can  justify  the  comparison 
made  in  Figs.  9a  and  9b  where  we  presented  the  temperature-dependent 
experimental  and  simulated  luminescence  rise  time.  Figure  11  shows  the 
luminescence  spectra  collected  at  the  band  edge  energy  region  when  the  electric 
field  of  the  light  is  0.6  x  10^  Vm“^  and  the  excess  photon  energy  is  50  meV.  The 
results  are  just  for  the  band  edge  energy  because  the  higher  energies  are  less 
important — in  this  case,  the  excess  photon  energy  is  lower  than  the  polar  optical 
phonon  energy. 
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Fig.  11  Calculated  luminescence  spectra  of  GaN  at  different  temperatures  collected  at 
band  edge  energy  region.  The  excess  photon  energy  is  50  meV  and  the  electric  field  of  the 
external  light  pulse  is  0.6x10’  V  m“*. 

The  solid  black  line  shows  the  results  when  the  temperature  is  50  K,  the  dashed 
red  line  shows  the  results  when  the  temperature  is  150  K,  and  the  dot-dashed 
green  line  shows  the  results  when  the  temperature  is  300  K.  As  seen  in  the  figure, 
when  the  temperature  decreases,  the  rise  time  of  the  luminescence  spectrum 
increases.  Figure  12  shows  the  results  of  the  integrated  luminescence  spectrum  in 
GaN  collected  up  to  8  ps  at  the  band  edge  energy  region.  The  lattice  temperature 
is  150  K,  and  the  excess  photon  energy  is  750  meV.  The  results  are  shown  versus 
the  normalized  power  of  the  laser  pulse.  The  solid  line  shows  the  results  of  the 
theory,  and  the  marker  points  show  those  of  experiment.  The  figure  clearly  shows 
a  reasonable  agreement  in  the  results.  As  the  power  increases,  because  of  the 
generation  of  more  carriers,  the  intensity  of  the  luminescence  spectrum  increases. 


Fig,  12  The  integrated  luminescence  spectra  of  GaN  collected  up  to  8  ps  vs,  normalized 
laser  power  collected  at  band  edge  energy  region.  The  excess  photon  energy  is  750  meV  and 
the  lattice  temperature  is  150  K. 
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2.5  Summary  for  Gallium  Nitride  (GaN)  Experiments  and 
Modeling 

We  presented  a  theoretical  and  experimental  study  of  the  subpicosecond  kinetics 
of  photoexcited  carriers  in  GaN.  We  considered  a  bulk  GaN  model  with  an 
isotropic  parabolic  heavy  hole  band  and  an  isotropic  nonparabolic  conduction 
band.  In  the  theoretical  model,  interaction  with  an  external  ultrafast  laser  pulse  is 
treated  coherently,  and  to  include  the  scattering  mechanisms  and  dephasing 
processes,  a  generalized  Monte  Carlo  simulation  is  used. 

The  scattering  mechanisms  included  are  carrier  interactions  with  polar  optical 
phonons  and  acoustic  phonons,  and  carrier-carrier  Coulomb  interactions.  The 
results  of  the  integrated  density  of  electrons  are  shown  for  2  different  excess 
photon  energies.  In  the  case  where  the  excess  photon  energy  (50  meV)  is  less  than 
the  polar  optical  phonon  energy,  carrier-carrier  scattering  is  the  dominant  process. 
However,  when  the  excess  photon  energy  (750  meV)  is  greater  than  optical 
phonon  energy,  the  polar  optical  phonon  scattering  is  the  most  significant  relaxing 
process.  Besides  the  carrier  densities,  temperature-dependent  luminescence 
spectra  have  been  studied.  We  showed  that  for  all  values  of  the  excess  photon 
energies,  as  the  temperature  decreases,  the  luminescence  rise  time  increases.  This 
effect  is  due  to  the  lower  polar  optical  scattering  rate,  which  results  in  the  fact  that 
this  process  is  no  longer  strong  enough  to  relax  the  carriers.  For  comparison,  we 
also  report  the  experimental  time-resolved  photoluminescence  studies  on  GaN 
samples.  In  conclusion,  we  have  shown  that  we  can  explain  the  kinetics  of  the 
photoexcited  carriers  in  GaN  by  including  only  carrier-carrier  and  carrier-phonon 
interactions  and  a  relatively  simple  2-band  electronic  structure  model. 
Furthermore,  we  have  presented  a  detailed  analysis  of  the  effect  of  a  different 
dephasing  machanism  on  the  carrier  kinetics.  Finally,  we  have  experimentally 
measured  and  simulated  using  our  theoretical  model  the  temperature  and  pulse 
power  dependence  of  the  photoluminescence  spectra. 

2.6  Results  and  Discussion  for  the  Zinc  Oxide  (ZnO)  Case 

In  this  section,  we  present  and  compre  the  normalized  luminescence  intensity 
results  obtained  by  theory  and  experiment  in  the  case  of  ZnO.  More  specifically, 
we  compare  the  luminescence  rise  times  and  the  effect  of  temperature  on  them. 
The  results  are  obtained  in  a  range  of  excitation  energies  and  temperature  values. 
The  material  parameters  for  ZnO  are  given  in  Table  2. 

The  calculated  polar  optical  scattering  rates  for  electrons  are  shown  in  Fig.  13, 
when  the  electric  field  of  the  light  is  0.6  x  10^  Vm  \  the  excess  photon  energy  is 
380  meV,  and  the  lattice  temperature  is  300  K. 


23 


Table  2  Parameters  used  in  the  simulation 


Parameter 

Value 

Unit 

771,. 

0.21772  0 

kg 

mi, 

0.59772  0 

kg 

(x 

•1.07f0 

Fni  ^ 

fs 

8.27ro 

Fm-i 

ho>Q 

72.9 

nieV 

E, 

3.373 

eV 

0 

M 

2 

A 

a 

0.4 

eV  ‘ 

D, 

7.5 

eV 

Dh 

22 

eV 

P 

6.10 

-3 

g  cm 

Vl 

6.59 

Ur  cm 

Vt 

2.79 

10*^  cm 
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cr 
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Fig.  13  Calculated  polar  optical  scattering  rate  for  the  nonparabolic  conduction  band  in 
ZnO  including  the  screening  effect.  The  solid  lines  show  the  phonon  emission  rate,  and  the 
dashed  lines  show  the  phonon  absorption  rate.  The  electric  field  of  the  light  is  0.6  x  10^  Vm“^, 
the  excess  photon  energy  is  380  meV,  and  the  lattice  temperature  is  300  K. 

Figures  14a  and  14b  show  the  results  of  the  normalized  luminescence  intensity 
obtained  by  experiment  and  simulation,  respectively.  The  excess  photon  energy  is 
380  meV,  and  the  lattice  temperature  is  300  K.  Different  curves  show  the  energy 
regions  where  the  luminescence  has  been  obtained,  and  each  curve  is  normalized 
to  its  maximum  value.  The  solid  black  line  shows  the  result  of  the  luminescence 
obtained  for  photon  energies  near  the  band  edge  transition.  The  dashed  red  line, 
dot-dashed  green  line,  and  double-dot-dashed  blue  line  show  the  same  results 
obtained  at  0.5,  1,  and  1.5  polar  optical  phonon  energy  in  excess  of  the  band  edge 
optical  energy  transition,  respectively. 
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Fig.  14  The  result  of  the  normalized  luminescence  intensity  from  ZnO  obtained  by 
experiment  (a)  and  simulation  (b).  The  excess  photon  energy  is  380  meV,  and  the  lattice 
temperature  is  300  K. 

Figure  15  shows  the  result  of  the  normalized  luminescence  intensity  from  ZnO 
obtained  at  the  band  edge  energy  region.  The  dashed  line  represents  the 
experimental  result,  and  the  simulated  result  is  shown  by  the  solid  line.  In  this 
case,  the  excess  photon  energy  is  40  meV,  and  the  lattice  temperature  is  75  K. 
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Fig.  15  The  result  of  the  normalized  luminescence  intensity  from  ZnO  obtained  by 
experiment  (dashed  line)  and  simulation  (solid  line).  The  excess  photon  energy  is  40  meV,  the 
electric  field  of  the  laser  pulse  is  0.6  x  10’  Vm“^,  and  the  lattice  temperature  is  75  K. 

Besides  the  results  of  the  rise  time,  we  also  investigated  the  effect  of  temperature 
on  the  rise  time.  Figure  16  shows  the  results  of  the  normalized  luminescence 
intensity,  when  the  excess  photon  energy  is  380  meV,  calculated  at  the  band  edge 
energy  region  over  a  range  of  lattice  temperature  values  from  150  to  300  K.  It 
clearly  shows  that,  as  the  lattice  temperature  increases,  the  luminescence  rise  time 
decreases.  This  is  because  at  lower  temperatures,  the  polar  optical  scattering  rate 
is  weaker  than  at  room  temperature  and,  therefore,  this  process  is  not  as  effective 
in  relaxing  the  carrier  distributions. 


Fig.  16  The  normalized  luminescence  intensity  from  ZnO  at  different  temperature  values 
calculated  at  band  edge  energy  region.  The  excess  photon  energy  is  380  meV,  and  the  electric 
field  of  the  external  light  pulse  is  0.6  x  10’  V  m“^. 
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In  this  study  we  considered  bulk  ZnO  modeled  by  an  isotropic  parabolic  heavy 
hole  band  and  an  isotropic  nonparabolic  conduction  band.  In  the  theoretical 
model,  interaction  with  an  external  ultrafast  laser  pulse  is  treated  coherently,  and 
the  scattering  mechanisms  and  the  dephasing  processes  are  taken  into  account 
through  a  generalized  Monte  Carlo  simulation.  The  scattering  mechanisms 
included  are  carrier  interactions  with  polar  optical  phonons  and  acoustic  phonons, 
and  carrier-carrier  Coulomb  interactions.  The  results  of  the  calculated  normalized 
luminescence  intensity  from  ZnO  at  2  excitation  energies  and  lattice  temperature 
values  are  in  good  agreement  with  those  of  the  experimental  study  obtained  by  the 
time-resolved  photoluminescence  method.  We  also  presented  the  calculated 
results  of  the  temperature-dependent  luminescence  intensity.  As  a  conclusion,  as 
the  temperature  decreases,  the  luminescence  rise  time  increases.  This  effect  is  due 
to  the  lower  polar  optical  scattering  rate,  which  results  because  this  process  is  no 
longer  strong  enough  to  relax  the  carriers. 

In  conclusion,  we  have  shown  that  we  can  explain  the  dynamics  of  the 
photoexcited  carriers  in  ZnO  by  including  only  carrier-carrier  and  carrier-phonon 
interactions  and  a  relatively  simple  2-band  electronic  structure  model. 

3.  Geometric  Single-Qubit  Quantum  Gates  for  an  Electron  Spin 
in  a  Quantum  Dot 


We  propose  a  scheme  to  perform  arbitrary  unitary  operations  on  a  single  electron- 
spin  qubit  in  a  quantum  dot.  The  design  is  solely  based  on  the  geometrical  phase 
that  the  qubit  state  acquires  after  a  cyclic  evolution  in  the  parameter  space.  The 
scheme  uses  ultrafast  linearly  chirped  pulses,  providing  adiabatic  excitation  of  the 
qubit  states,  and  the  geometric  phase  is  fully  controlled  by  the  relative  phase 
between  pulses.  The  analytic  expression  of  the  evolution  operator  for  the  electron 
spin  in  a  quantum  dot,  which  provides  a  clear  geometrical  interpretation  of  the 
qubit  dynamics,  is  obtained.  Using  parameters  of  InGaN/GaN  and  GaN/AlN 
quantum  dots,  we  provide  an  estimate  for  the  time  scale  of  the  qubit  rotations  and 
parameters  of  the  external  fields. 

3.1  Introduction 

The  electron  spin  in  a  single  quantum  dot  is  one  of  the  perspective  realizations  of 
a  qubit  for  the  implementation  of  a  quantum  computer.  During  the  last  decade, 
several  control  schemes  to  perform  single  gate  operations  on  a  single  quantum  dot 
spin  have  been  reported.  Here  we  propose  a  scheme  that  allows  performing 
ultrafast  arbitrary  unitary  operations  on  a  single  qubit  represented  by  the  electron 
spin.  The  idea  of  geometric  manipulation  of  the  qubit  wave  function  has  been 
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recently  developed  into  a  new  research  direction  called  geometric  quantum 
computing.^^"*^  The  main  motivation  of  this  development  is  the  robustness  of 
geometric  quantum  gates  against  noise. In  this  report,  we  demonstrate  how  to 
use  the  geometric  phase,  which  the  Bloch  vector  gains  along  the  cyclic  path,  to 
prepare  an  arbitrary  state  of  a  single  qubit.  We  show  that  the  geometrical  phase  is 
fully  controllable  by  the  relative  phase  between  the  external  fields.  We  also 
discuss  the  realistic  implementation  of  the  proposed  design  using  the  electron  spin 
in  a  charged  quantum  dot  as  an  example  of  a  qubit. 

3.2  General  Equations  of  Motion 

Let  us  consider  the  coherent  Raman  excitation  in  the  3-level  A -type  system 
consisting  of  the  2  lowest  states  of  electron  spin  1 0)  M  -X)  and  1 1)  MX)  coupled 
through  an  intermediate  trion  state  IT)  consisting  of  2  electrons  and  a  heavy 
hole^i  (Fig.  17). 

'i  "■■■V' . 

-t-L|i)  =  |x) 

-fd— |o>^|-v> 

Fig.  17  Energy  structure  of  the  3-level  system  comprising  the  2  electron  spin  states  and  the 
trion  state 

We  assume  that  the  trion  state  is  far  off-resonance  with  the  external  fields  to 
ensure  that  decoherence  on  the  trion-qubit  transitions  can  be  neglected.  The 
electron  spin  states  are  split  by  an  external  magnetic  field;  the  separation  energy  is 
ho)^ .  The  total  wave  function  of  the  system 

I  'P(O)  =  a^{t)  1 0)  +  1  +  I T),  (17) 

where  aoj(t)  and  are  the  probability  amplitudes,  is  governed  by  the  time- 
dependent  Schrodinger  equation  with  the  Hamiltonian 

^  0  0  -(Q^(0  +  Os(0)^ 

n  =  h  0  0}^  -(Qpio+nsiO)  ,  (18) 

^-(Qp(r)-i-Q^(r))  +  coj  ^ 

where 
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5  ^po,so  cos[(yp  5 

^^P,S  ^P0,50  ^F,5 

^po,so^^^  ^0T,Ti^p,s  ^  ^  5  3nd 

Opo,so(^)  =  Mit,to^p,s  (^)  /  ^  the  Rabi  frequencies; 


/^07-,n  are  the  dipole  moments; 

^'p  5  (t)  are  the  pulse  envelopes; 

^p.s  are  the  center  frequencies; 

^z^p  ^CO  are  the  time-dependent  phases;  and 
TicOj  is  the  energy  of  the  trion  state. 


We  are  considering  here  a  case  of  linearly  chirped  pulses  such  that 


^P,S  “  ^P,S  ^ 

where  <f)p  ^  are  the  initial  phases  and  ap  ^  are  the  chirps  of  the  pulses. 
Using  transformation  I  'P(t))  =  I  ,  where 


U 


RWA 


^1  0 

Q  ^i(cos-o)p)t 


0 


= 

RWA 


0 


0 


0  ^ 
0 

-iCDp. 

0  ^ 


-1  —I  0  Q 

^0  0 


(19) 


(20) 


(21) 


we  make  the  rotating  wave  approximation  (RWA)  by  neglecting  the  rapidly 
oscillating  terms  with  frequency  2(Dg,2(Op  dxvdcOg+cOp.  In  the  RWA,  the 

Hamiltonian  has  the  following  form: 
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H  =U-UHU 


RWA  RWA  Usm 


+n 


h 

2 

0 

0 

0 

V 

h 

2 


0 

0 


Qp_^e 


0 

-2®, 


0 

-Ao) 


0  ^ 
0 


0  -(e) 

0 
0 


pj 


Q*  I  r^*  J^cot 


0 

2(A(i)-(i)J 


Qp,e 


Qp^ay'^^'+Q 


-iAcut  \ 


S+ 


S+ 


-2(0t 


n.. +«*.«■“"  ' 


-5  + 


5  + 


-2  A, 


(22) 


where 

Ap=(Oj  —(Op,  A(o  =  (Op  —(Og , 


Qp,=Qpo(Oe'^^^^Q,,=Q,oay 

Qp^=Qpo(Oe'^^'\and 


¥s(0 


Qc+  —  ^co(^)^ 


'Vx  (0 


Assuming  large  detunings  of  the  pump  and  Stokes  field  frequencies  from  the 
transition  frequencies  to  the  trion  state,  we  apply  the  adiabatic  elimination  of  the 
trion  state.  In  that  case,  ^  «  0  and  we  find 

b  =  +  (23) 

Substituting  Eq.  23  into  equations  for  Oq  j  after  some  algebra,  we  obtain  the 
following  form  of  the  Hamiltonian  for  the  effective  2-level  system: 


H 


(24) 


where  S  =  (o-  Aco , 


^aco(^)  “  .  (^po(^)  +  ^L(^)  +  2Qpg(t)Q^g(t) COS \(f>p{t)  —  (f)^{t)  +  A<»t],  (25) 

ZjL\p 

^aci  (0  =  :r^  (f^L  (0  +  (0  +  20^0  (O^po  (0  cos  [(/)p  (0  -  4  (0  -I-  A(0t\,  (26) 

Z^L\  n 
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2A, 


(27) 


+  Q^o  (OO50  (Oe''^"  +  Opo  (O^^so  )• 


Since  QpQ(f)Q5Q(/)  =  Qpq(?)Q^q(?)  ,  we  have 

-iAcot 


(28) 


+  2Qpo  (?  )Q5o  (0  cos  [^zip  (? )  -  4  (0  +  A<»?] ). 
Making  the  transformation  I  'P(t))  =  U I  'P(t)) ,  with 


U  = 


-idtl2+i[  dtQ.^f(t) 
o  J-00 


0 


-idtl2+i[  dtQ.^f(t) 

n  J-00 


-iSt/2+i  [  dt  (r  ) 

=  e  •'-”  I, 


(29) 


where  Q^j(r)  =  (f2^co(^)  +  ^aci(^))(4,  we  ean  rewrite  Eq.  24  in  a  more  symmetric 
form.  In  the  new  basis,  the  Hamiltonian  takes  the  form 


H  =  U"‘HU  -  U  =  H  +  ^(Q^,  (0  -  ^  /  2)  I 


eff' 


(30) 


where  Q^,y  (0  =  (Q,,o(0  -  ^acAt))  >  2 . 


Taking  into  aceount  the  exact  time  dependenee  of  the  phases  in  Eq.  19,  we  apply 
another  transformation  I  'F(r))  =  I  'E(r))  ,  with 


0 


0 


J 


(31) 


where  <^  =  .  Therefore,  in  the  field  interaetion  representation  the 

Hamiltonian  takes  the  following  form: 
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-5-n,,^{t)y  -a 

2[  -5+0-^,, it)^ 

^_nf_  Sit) 

2[n*^(t)e-‘^^  -Sit)  / 

where  Sit)  =  S  +  ia^  -ap)t  +  it) ,  S<f>  =  ^p-^^,  and 


^eff  (0  - 


Q^oCOQ^oCpr  ^  Q^o(OQ^o(0  +  Q,o(OQ,o(0 

^-i{lt^(l)+ls.cot+^f'  H 

2^p  y  Qpo(OQso(0 

I  ^PO  (OQ^o(0 

^-2i{^(|)+^cot+<^t^  12^ 

flj,o(0^so(0  j 


(33) 


In  general,  the  differential  AC  Stark  shift,  it) ,  is  not  zero  and  has  to  be  taken 

into  account.  Using  the  definition  of  the  Rabi  frequency,  we  can  rewrite  Eqs.  25 
and  26  in  the  form 

2 

^acoiO  =  Ay-  (Elit)  +  Elit)  +  2Epit)E,it)  cos  +  Sat],  (34) 

2n  Ap 

2 

^«ci(0  =  A\  {Elit)  +  El  it)  +  2Epit)E^  it)  cos  [(l)pit)  -  4(0  +  Scot],  (35) 

In  Ap 

In  the  case  when  «  ju^j ,  the  AC  Stark  shifts  are  the  same  for  both  qubit  states, 
so  that  (?)  =  0 ,  and  we  have 


2\vi]ffit)e-^^^  -S  +  Ct  ; 


(36) 


Note  that  we  do  not  require  completely  overlapped  pulses  here. 


3.3  Nonimpulsive  Case 

In  some  excitation  schemes,  because  of  selection  rules  that  take  into  consideration 
the  polarization  of  the  external  field,  the  pump  and  Stokes  field  interact  only  with 
the  corresponding  transitions,  and  the  general  Hamiltonian  can  be  simplified.  For 
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that  case,  we  can  obtain  the  correct  form  of  the  Hamiltonian  by  putting 
Qpo(t)  =  ^so(0  =  0  in  Eqs.  25,  26,  and  33.  It  results  in 


^dif  (0 


^[Q^o(0-QL(0], 

4Ap 


^eff  (0  - 


^Po(0^so(0 

2Ap 


(37) 

(38) 

(39) 

(40) 


Therefore,  in  the  field  interaction  representation,  the  Hamiltonian  has  the  form 

hf  d(t) 

2  (Oe 


H 


-S(0 


(41) 


where 


<5(t)  —  <5  +  (oTc  —  (Xp)t  +  [f2pQ(t)  ~  f2jQ(r)]  /  4Ap. 


(42) 


3.4  Adiabatic  Solution 


The  Hamiltonian  in  Eq.  41  controls  the  dynamics  of  the  qubit  wave  function  in 
the  approximation  of  the  adiabatic  elimination  of  the  trion  state.  Here  we  consider 
the  adiabatic  excitation  of  the  qubit  and  find  the  adiabatic  solution  of  the 
Schrodinger  equation  with  the  Hamiltonian  in  Eq.  41. 

Since  the  phase  factor,  ,  of  the  coupling  term  in  Eq.  41  is  time  independent,  it 
is  convenient  to  use  the  following  transformation  1 0(t))  =  A I  'F(O)  ,  where 

A  =1 0)<0 1  +e'^'^  1 1><1 1=  ,  (43) 


so  that  the  new  wave  function  is  governed  by  the  Hamiltonian 


H=AHA  = 


h 


5{t)  0,(0' 

vO,(0  -^(0, 


h 


((7(0  a, +0,(0  oj.  (44) 


To  solve  the  Schrodinger  equation  in  the  adiabatic  representation,  we  apply 
another  transformation,  1 0(r))  =  R(r)  1 0(r)) ,  where 


R(0  = 


cos6>(0  sin6>(0 

-sin^(0  cos  6^(0 


=  e 


i6(t)a 


(45) 
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and  tan[2^(?)]  =  Qg(?)/  S{t) .  In  the  new  basis  the  Hamiltonian,  Eq.  44,  takes  the 
form 


H(0  =  R  WH(f)  R^‘d)  = 


0 


0  1 


-  -hA  (t)  0^/2  , 


(46) 


where  2  (t)  =  +  . 


As  we  see,  the  Hamiltonian  in  Eq.  46  is  diagonal  in  the  adiabatic  basis,  and  we 
can  readily  write  down  the  solution.  However,  since  the  transformation  R(Ois 

time  dependent,  an  additional  nonadiabatic  coupling  term  is  present  in  the  general 
Schrodinger  equation 

=  a\0{t))-hd{t)a\^(t)),  (47) 

2 


where 


2(Ql(t)  +  S\t)) 


(48) 


Neglecting  the  nonadiabatic  coupling  term  in  Eq.  47,  we  readily  obtain  for  the 
qubit  wave  function  in  the  original  basis 


I  ^(0)  =  e  2  ' 


e  ^ 


W(0)ay 

e  ^ 


e  ^  \  ^(0)>  , 


(49) 


where  A(0  =  dt'A(t')  .  The  general  form  of  the  evolution  operator  in  Eq.  49  is 

Jo 

well  justified  if  the  following  condition  IQ^(0<^(0“^^e(0^(0D  ^^(t)  is  valid. 


In  the  case  of  completely  overlapped  pulses,  ,  with  identical  chirp 

rates,  ap=a^,  for  the  resonant  qubit,  S(t)  =  S  =  0,  we  have  9(t)  =  9(0)  =  7rl 4  , 
and  the  transformation  matrix  becomes  R(r)  =  R(0)  =  .  Therefore,  the 

unitary  evolution  operator  for  the  wave  function  of  the  resonant  qubit  takes  the 
form 

J  cos  (5(0/ 2)  sin  (5(0/ 2)^ 

^*'^^~[ie-‘^^sin(5(0/2)  cos (5(0/ 2)  ^ 

=  cos  (5(0  /  2)  I  +  /  sin  (5(0  /  2)  (e‘^^a^  + 

=  cos  (5(0  /  2)  I  -  /  sin  (5(0  /  2)  (n-o)  =  ^  (50) 
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where  5(0=1'  )  effective  pulse  area,  n  =  (-cosA^,sinA^,0) , 

o+  =  (o^  ±/o^)/2are  the  Pauli  raising  and  lowering  operators,  and  are  the 

Pauli  operators.  The  nonadiabatic  coupling  term,  Eq.  48,  is  zero  for  the  resonant 
qubit,  and  the  solution  of  the  Schrodinger  equation  in  the  adiabatic 
approximation,  Eq.  50,  is  the  exact  solution. 


The  density  plots  of  the  population  and  coherence  ( I  aQ(7’)a*(7’)  I )  at  the  final  time 
(after  the  pulse  excitation)  as  a  function  of  the  effective  pulse  area,  S(T) ,  and  the 
dimensionless  frequency  chirp  parameter,  a  /  tI,  described  by  the  unitary 
evolution  operator  in  Eq.  50  is  depicted  in  Eig.  18.  We  use  the  Gaussian  shape  for 
the  pulse  envelopes,  assuming  that  linear  chirp  is  obtained  by  applying  a  linear 
optics  technique,  meaning  that  a  transform-limited  pulse  of  duration  Tq  is 

chirped,  conserving  the  energy  of  the  pulse.  The  temporal  (a)  and  spectral  a  ^ 

chirps  are  related  as  a  =  aT'^  ! (\  +  a^ ! ,  where  is  the  transform-limited 

pulse  duration.  We  observe  the  Rabi  oscillation  regime,  when  the  population  of 
the  qubit  states  is  changing  between  0  and  1  while  the  coherence  is  changing 
between  0  and  1/2.  This  behavior  does  not  depend  on  the  chirp  rate,  since  the 
effective  Rabi  frequency  is  determined  by  the  product  of  the  pump  and 

Stokes  Rabi  frequencies:  Qpo5Q(0  =  nQexp{-r^/(2r^)}/[l-l-ci'^/ro]'^'‘with  the 
chirp-dependent  pulse  duration  r  =  Tq\\  +  and  the  amplitude.^^’^^ 


Fig.  18  The  density  plot  of  the  1 1)  state  population  (a)  and  coherence  (b)  as  a  function  of 
the  effective  pulse  area  and  frequency  chirp;  ,  s  =  o.  Initially,  only  the  1 1)  state  is 

populated. 

In  turn,  for  the  off-resonant  qubits,  S  the  evolution  operator  in  the  adiabatic 
approximation  takes  the  form 
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U(0  = 


(51) 


^  cos  0{t)  sin  ^(0^ 

sin  e{t)  cos  0(t)  ^ 


where 


cos^(0  =  ^=  |1  + 


S 


^/2^  ^s^+nl(t) 


,  sin  0(t)  = 


^/2  V  p^+nl(t) 


(52) 


and  <^(0  =  -  £  p^+Ql(t')dt' 


is  the  effective  pulse  area. 


For  the  off-resonant  case,  the  transformation  matrixes  are 

^cos6>(0  -sin6>(t)^  _  _  f  1  0^ 

^sin6'(t)  cos6>(0  y 


R-'(0- 


R  (0): 


0  1 


(53) 


Figure  19  demonstrates  excitation  of  the  off-resonant  qubit.  As  expected,  the 
population  of  the  off-resonant  qubit  at  the  final  time  is  not  changed  by  the 
external  fields  as  long  as  the  pulse  exeitation  parameters  are  in  the  adiabatie 
regime.  This  is  the  regime  of  adiabatic  return.  However,  we  still  observe  the  Rabi 
oscillation  for  the  value  of  the  chirp  I  a  l<  5rQ .  This  is  the  area  of  nonadiabatic 

population  transfer  where  the  nonadiabatic  coupling  term  cannot  be  neglected  in 
Eq.  47. 
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Fig.  19  The  density  plot  of  the  1 1)  state  population  (a)  and  coherence  (b)  as  a  function  of 
the  effective  pulse  area  and  frequency  chirp;  St^  =  0.75 .  Initially,  only  the  1 1)  state 

is  populated. 
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3.5  Bloch  Vector  Representation 


The  dynamics  of  the  qubit  wave  function  can  be  described  equally  well  using  the 
Bloch  vector  representation.  In  addition,  the  Bloch  vector  formalism  allows  a  very 
nice  and  clear  geometrical  interpretation  of  qubit  dynamics.^'^  In  this  section,  we 
give  a  short  overview  of  the  Bloch  picture. 

A  general  state  of  a  qubit  can  be  described  as 


I  'F>  =  ^0  (0 1 0>  +  (0 1 1>  =  cos  (y5  /  2)  1 0>  +  sin  (y5  /  2)  1 1>,  (54) 


where  a  and  are  the  phase  parameters.  Up  to  an  insignificant  global  phase, 
the  wave  function  can  be  mapped  into  a  unitary  Bloch  vector  B  =  (m,  v,  w)  ,  as 
shown  in  Fig.  20. 


t  z 


X 


Fig.  20  The  Bloch  vector  representation  of  the  qubit  state.  Excitation  of  the  qubit  by  an 
external  field  corresponds  to  the  rotation  of  the  B  vector  about  the  pseudo  field  vector,  O, 
with  components  determined  by  the  effective  Rabi  frequency  detuning  5  ^  and  the 

relative  phase  A^. 

To  use  Bloch  vector  representation,  we  construct  the  qubit  density  matrix 


vflo  Qn)  '2.\e'“smj3  l-cos/?J’ 


(55) 


where  =  a.  (t)a*j  (t) ,  i,  j  =  0,1.  Using  a  Pauli  matrix  decomposition 


(56) 
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we  identify  relation  between  the  components  of  the  qubit  wave  function,  the  qubit 
density  matrix  elements,  and  the  Bloch  vector  components.  Therefore,  we  obtain 
the  following  expression: 


Qoi  ^10 

^cosasin 

B  = 

V 

= 

KQoi  ~  Qio) 

= 

sin  a  sin  /? 

^  Qoo~Qii  J 

^  cos  ^ 

Taking  into  account  the  Hamiltonian  in  Eq.  41,  the  equation  of  motion  for  the 
density  matrix  /^p  =  H  p— p  H  takes  the  following  explicit  form 

'fti =-iS(')ai  -(01  -ao)f2,  (<)«“'  /2, 

(5o) 

100  =  <^(1)00  +(01  -  /  2, 
;a,=-(e-'«ai-e'«a„)n,<0/2. 

Using  the  relations  in  Eq.  57,  we  can  also  write  the  dynamic  equation  in  the  Bloch 
vector  representation  as 

u  =  vS(t)  +  wQg(0sin  A^, 

v  =  -uS(t)  +  wQ.^(t)cosA(f>,  (59) 

w  =  -uD.^  (t)  sin  (t)  cos  A^. 


Introducing  a  pseudo  field  vector,  Q,  with  components  determined  by  the 
effective  Rabi  frequency,  2-photon  detuning,  and  the  relative  phase  between 
pump  and  Stokes  pulses. 


a  = 


^-Q^(0cos  A^^ 
Q/OsinA^zi  , 


I  -^(0 


J 


(60) 


we  can  rewrite  Eq.  59  in  the  compact  form 

B  =  nxB.  (61) 


This  is  the  Bloch  equation,  which  describes  a  precession  of  the  Bloch  vector,  B, 
about  the  pseudo  field  vector,  Q,  and  allows  clear,  intuitive  interpretation  of  qubit 
dynamics. 
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3.6  Evolution  Operator  of  the  Bloch  Vector 


Since  we  already  know  the  exact  form  of  the  evolution  operator  U(t)  in  Eq.  50 
and 


l'P(0)  =  U(0l'P(0)), 


(62) 


we  can  easily  construct  the  evolution  operator  for  the  Bloch  vector.  Using  the 
definition  of  the  density  matrix,  we  have 

P(0  =1  'P(0><'P(0 1=  U(0 1  'P(0)><'P(0)  I  U^(0  =  U(0  p(0)  r  (0,  (63) 


where 


p(0)H'E(0)><'F(0)h 


^*00  ((^)  031  ((^) 

,0o((^)  0i(O), 


1  +  Wo  u^-iv^ 
\Ho  +  i^o  l-%  y 


is  the  initial  condition  and 


t^A  - 


u^(0 


cos(5(0/2)  -ie'^'^sm{S(t)/2) 

cos  (5(0/ 2)  ^ 

=  cos  (5(t)  /  2)  I  -  /  sin  (5(0  /  2)  (e'^o^  +  V) 

=  cos  (5(0  /  2)  I  +  /  sin  (5(0  /  2)  (n-o)  =  . 


Therefore, 


p(0  =  u(0p(0)U00 


trA  - 


cos  (5(0/ 2)  sin  (5(0/ 2)' 

sin  (5(0/ 2)  cos  (5(0/ 2) 


^0oo(O)  ^oi(0)Y  cos(5(0/2)  -^V"^in(5(0/2)^ 


,0o((^)  0ii(O), 


-ie-‘^^sm{S(t)/2)  cos  (5(0/ 2) 


Taking  into  account  Eq.  57,  we  obtain  for  the  Bloeh  vector: 

^C^  +  S^cos{2A(p)  -5"sin(2A^^)  205  sin  (A^^)^ 

B(0  =  |  -5"sin(2A<zi)  C"  -  5"  cos(2A<zi)  2C-5  cos  ( A^zi) 

-2G5sin(A^^)  -2C-5cos(A<;i)  C"-5" 


(64) 


(65) 


B(0),  (66) 


where  C  =  cos  (5(0/ 2),  5  =  sin  (5(0/ 2). 


3.7  Ultrafast  Qubit  Rotations  Using  Geometrical  Phase 

At  this  point,  we  are  ready  to  discuss  implementation  of  the  single  qubit  gates 
sinee  we  have  obtained  the  anal3^ie  solution  for  the  qubit  wave  function  and 
constructed  the  evolution  operator  in  the  Bloch  vector  representation.  A  universal 
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set  of  quantum  gates  has  been  intensively  discussed  in  the  literature  related  to  the 
universality  in  quantum  computation.^^  To  perform  quantum  computation,  we 
must  have  2  major  building  blocks  at  our  disposal:  arbitrary  unitary  operations  on 
a  single  qubit  and  a  controlled-NOT  operation  on  2  qubits.  Here  we  address  only 
single  qubit  manipulation. 

To  demonstrate  arbitrary  geometric  operations  on  a  single  qubit,  we  use  the  Bloch 
vector  representation  discussed  in  the  previous  section.  Since  any  unitary  rotation 
of  the  Bloch  vector  can  be  decomposed  as 

R^(a,)  R/a,)  R^(a,),  (67) 


where  R,.  =  e'""'  (i  =  y,z)  are  the  rotation  operators,  we  need  to  demonstrate 

rotations  of  the  qubit  Bloch  vector  about  the  z  and  y  axes  by  applying  various 
sequences  of  external  pulses.  The  decomposition  in  Eq.  67  plays  an  important  role 
in  circuit-based  quantum  computing,  as  it  shows  explicitly  that  2  single-qubit 
operations  allow  us  to  create  the  arbitrary  state  of  the  qubit.  Here  we  show  how 
this  can  be  accomplished  by  controlling  the  parameters  of  the  external  pulses, 
which  are  defined  by  the  explicit  form  of  the  evolution  operator  (see  Eq.  50). 
There  are  2  distinct  ways  of  the  implementation  of  rotation,  depending  on  which 
part  of  the  total  qubit  phase  we  employ:  dynamical  or  geometrical.  Quantum  gates 
relying  on  geometrical  quantum  phases  are  called  holonomic  gates,  and  they  are 
expected  to  be  robust  with  respect  to  noise.^^  *^ 


To  implement  the  rotation  of  the  Bloch  vector  about  the  z  axis  (the  phase  gate) 
based  on  the  geometrical  phase,  we  can  use  the  evolution  operator  of  the  resonant 
qubit,  Eq.  50.  The  product  of  2  evolution  operators  corresponding  to  the  sequence 
of  2  pulses  with  the  relative  phase  =  ^ + ;r  gives 


0 


0 


-l(p 


(68) 


where  the  first  subindex  of  U  indicates  the  pulse  area,  S(T) ,  and  the  second  one 
indicates  the  phase,  . 


Figure  21  shows  the  Bloch  vector  trajectories  of  the  qubit  basis  states  1 0)  and  1 1), 
which  correspond  to  the  angles  j^  =  0  and  P  =  7l  in  Eq.  57  and  the  Bloch  vector 
initially  pointing  in  z  and  — z  directions  while  the  vector  nj=(-Q^,0,0)  is 
pointing  in  -jc  direction.  For  simplicity,  we  chose  A^  =  0  for  the  first  n:  -pulse. 
The  first  ;r -pulse  flips  the  population  to  the  state  II)  (10));  correspondingly,  the 
Bloch  vector  turns  about  the  effective  field  vector  (about  the  x  axis),  and  it 
stays  in  the  y,  z  plane  all  the  time  and  points  in  the  — z  ( z )  direction  at  the  end  of 


40 


the  pulse.  Because  of  the  second  ;r -pulse,  the  population  is  transferred  back  to  the 
initial  state  10)  (II));  therefore,  the  Bloch  vector  returns  to  its  original  position 
pointing  along  the  z  (— z)  axis.  However,  since  we  chose  b,^  =  (p-\-7t  for  the 

second  n  -pulse,  the  pseudo-field  vector  is  rotated  counterclockwise  by  the  angle 
in  the  x,  y  plane,  Qj  =  0)  ’  the  Bloch  vector  moves 

in  the  plane  perpendicular  to  the x:,y plane  and  has  the  angle  TiH-a 
{-7t 1 2- a)  with  the x, z  plane. 


Fig.  21  The  Bloch  vector  trajectory  for  the  quhit  state  1 0)  in  panel  (a)  and  the  quhit  state 
1 1)  in  panel  (h)  generated  hy  the  sequence  of  2  ;'r  -pulses  with  the  relative  phase  (p+K 


The  Bloch  vectors  representing  a  pair  of  orthogonal  basis  states  1 0)  and  1 1) 
follow  a  path  enclosing  correspondingly  solid  angles  of  2(p  and  -2(p .  The 

geometrical  phase  is  equal  to  one-half  of  the  solid  angle,  which  means  the  basis 
states  1 0)  and  1 1)  gain  phases  (p  and  -cp ,  and  the  evolution  operator  takes  the 

form  of  the  phase  gate,  Eq.  68,  with  the  relative  phase  controlling  the  phase  of  the 
gate. 

The  rotation  operator  about  the  y  axis  can  be  constructed  using  3  pulses.  The  first 
and  third  pulse  is  /  2  -pulse  with  =  0 ,  while  the  second  pulse  is  k  -pulse  with 
the  relative  phase  n  +  cp.  It  is  easy  to  show,  using  Eq.  50,  that  this  3-pulse 
sequence  results  in 


R  (^)  =  U,  U, 


-;0 

2 


n\n^(p 


-;0 

2 


(69) 


To  demonstrate  the  geometrical  nature  of  the  R^(^3)  operation,  we  use  the  fact 

that  it  creates  the  relative  phase  between  the  qubit  basis  states 
I  +/)  =  (1 0)  ±  / 1 1))  /  a/2  .  In  the  Bloch  representation,  these  states  have  the  form 


±i- 

COS 

UJ 

1 0)  -1-  e  ^  sin 

UJ 

II), 
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which  are  2  vectors  defined  by  the  angles  f5  =  7t  12  and  a  =  ±n !  2  and  pointing 
in  the  y  and  -y  directions,  as  shown  in  Fig.  22.  The  trajectory  of  the  Bloch 
vector  representing  the  states  I  +/)  is  shown  in  Fig.  22.  The  pseudo-field  vectors 
and  are  defined  by  the  effective  Rabi  frequencies  of  the  first  and  third 
pulses  and  are  pointing  in  the  -jc  direction  since  =  0 .  The  second  pseudo¬ 
field  vector  is  rotated  counterclockwise  by  the  angle  in  the  x,y  plane, 

the  same  as  in  the  case  above.  The  initial  Bloch  vector  is  pointing  in  the  y  (-y) 
direction.  The  first  tt/ 2  -pulse  rotates  the  Bloch  vector  about  flj  to  the  position 
of  the  state  II)  (10)).  The  second  pulse  flips  the  direction  of  the  Bloch  vector. 
The  third  ;t/ 2 -pulse  returns  the  Bloch  vector  to  its  original  position.  The  Bloch 
vector  and  the  pseudo-field  vector  are  orthogonal  during  the  whole  evolution. 
Similar  to  the  previous  case,  we  observe  that  the  basis  states  I  i)  and  I  —i)  follow  a 
path  enclosing  correspondingly  solid  angles  of  2(p  and  -2(p .  Therefore,  they  gain 
the  relative  phase  2(p,  which  is  the  geometrical  phase  defined  by  the  relative 
phase  between  pulses.  It  is  easy  to  show  that  the  phase  gate  in  the  I  +/)  basis  is 
equivalent  to  the  Ry(^£>)  gate  in  the  10),  II)  basis. 


Fig.  22  The  Bloch  vector  trajectory  for  the  quhit  state  I  i)  in  panel  (a)  and  the  qubit  state 
I-/)  in  panel  (b)  generated  by  the  sequence  of  2  ;r/ 2 -pulses  and  one  Tt  pulse  with  the 
relative  phase  (p+Ti 


3.8  Rotation  in  the  Bloch  Representation 

Rotation  operations  in  the  Schrodinger  picture  are 

l'i'>  =  R,IT„>, 


where  i  =  {x,y,z} ,  I  'Fg)  is  the  initial  wave  function. 


R 


^  cos  <p  i  sin  (p^ 
J  sin  <p  cos  <p  ) 


(70) 
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R 


_  ”)• 


^  cos^  sin^^ 
^-sin^  CO?,  (pj 


R 


r^iv 

0 


0 


y 


Therefore,  for  the  density  matrix,  we  have 

pH'FX'FNR,  I'FoX^JR:^  =R,  k:\ 

where  Pg  =1  'Pq)('Fq  I  is  the  initial  density  matrix. 


(71) 


Using  equation 

B  =  rr[o  p]  =  (rr[(T^  9\,Tr[G^  p],7’r[o^  p])  =  (M,v,w), 


we  find  the  following  expressions  for  the  transformation  of  the  Bloch  vector 
components 

M  =  rr[o^R.PoR-‘], 
v  =  rr[OyR,.PoR:‘], 
w  =  rr[c^R,.poR"']. 


We  obtain 


R 


R 


R, 


^10  0  ^ 
0  cos (2^)  sin (2^)  , 
^0  -sin  (2^)  cos  (2^)^ 

^cos(2^)  0  -sin(2^(?)^ 
0  10, 
^sin(2^)  0  cos(2^)^ 

^  cos(2^)  sin(2^)  0^ 

-sin (2^)  cos (2^)  0  . 


(72) 


(73) 


(74) 


The  evolution  operator  in  the  Bloch  representation 


U(0  = 


C"  +  5"cos(2A^)  -5"sin(2A^)  2G5sin(A^) 

(2A^?i)  C"-S"cos(2A<^)  2C-Scos(A<^) 


-S^  sin 


-20  S  sin  (  A^)  —20S  cos  (  A^) 


C  -S' 


(75) 


43 


Therefore,  to  demonstrate,  for  example,  the  rotation  in  the  Bloch  picture,  we 
see  that  the  first  n  pulse  with  the  relative  phase  =  0  gives  for  the  evolution 
operator 


U 


;t;0 


^10  0^ 
0-10 
lo  0  -ij 


(76) 


The  second  n  pulse  with  the  relative  phase  t^(j)  =  n (p  gives  for  the  evolution 
operator 


U 


7l\7l-\-(p 


^  COS {2{7r  +  (p))  -sin(2(;T  +  ^))  0^ 

-sin(2(;^  +  ^))  -cos(2(;r  +  ^))  0 


(77) 


Finally,  the  sequence  of  and  q  results  in 


''  cos(2(;t-i-^)) 

-sin(2(;T-l-^)) 

0^ 

fl 

0 

0^ 

7i:\7i:+(^  7i:\0 

-sin(2(;T-l-^)) 

-cos(2(;t-i-^)) 

0 

0 

-1 

0 

.  0 

0 

-ij 

vO 

0  -ij 

^  cos(2^)  sin  (2^)  0^ 

-sin  (2^)  cos(2^p)  0 


(78) 


Similarly,  one  could  consider  the  R  rotation  in  the  Bloch  picture. 


3.9  Generalization  of  the  Single-Qubit  Operation  Using  Bright- 
Dark  Basis 

In  the  previous  sections,  we  have  considered  several  excitation  schemes  of  the  3- 
level  system  and  discussed  a  possible  implementation  of  single -qubit  gates.  It  was 
shown  that  all  possible  qubit  states  can  be  created  in  a  controllable  fashion  using  a 
couple  of  completely  overlapped  laser  pulses,  Q.p(t)  =  Q.g(t) .  In  this  section,  we 

present  a  more  general  solution,  which  allows  some  additional  flexibility  in  terms 
of  the  ratio  of  the  pump  and  Stokes  pulse  amplitudes.  Again,  we  address  here  the 
coherent  Raman  excitation  in  a  3-level  A  -type  system  consisting  of  the  2  lowest 
states  of  electron  spin  1 0)  M  -X)  and  1 1)  M  X)  coupled  through  an  intermediate 
trion  state  I T)  (see  Fig.  17)  and  assume  that  the  trion  state  is  far  off-resonance 
with  the  external  fields.  In  addition,  we  restrict  our  consideration  to  the 


44 


nonimpulsive  regime  and  can  then  put  =  =  0  in  Eq.  22  so  that  the 

Hamiltonian  takes  the  form 


H  H  ih  UwA 


h 

^  0 

0 

"0 

0 

0  ^ 

2 

0 

~2<Ob 

+  fi 

0 

-Aft) 

0 

—2cOjj 

.0 

0 

—(Op  j 

^  0  0  " 
0  2(A&)-®p)  , 

q;,  -2Ap^ 


where 

t^p  =  (Oj  —  (Op  , 

Ao  =  (Op—(o^, 

Q.p,=Q.p,{ty*^^^\ 

Q,,=Q,o(ty^^«,and 

(f>p  g  (t)  =  (f>p  g  +  oCp  /  2 . 


(79) 


Let  us  consider  the  case  when  Q.pQ(t)  =  D.Q(t)  cos  &  and  Q^qCO  =  i2Q(0sini9  and 

keeping  the  same  time-dependent  envelope  for  the  pump  and  Stokes  Rabi 
frequencies  while  the  mixing  angle  3  controls  the  ration  between  the  maximum  of 
the  Rabi  frequencies. 


In  the  bright-dark  basis,  the  Hamiltonian  takes  the  form 


h=r„hr;;  =  * 


2SS^ 

15SCe'** 

^-^:(OQo(0 


15SCe'^^ 

2dC^ 


(80) 


where 

A  =  sin  <9  5 


C  =  cos  i9  5 

~  t 

K  />\  iat^  12/^2  .  iBt^  12  ri2  -• 

^^{t)  =  e  C  +e'^  S  ,and 
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In  the  resonant  case,  =  0 ,  we  obtain 

r  0  0  -^,(OQo(o 

H  =  R,,H  R;‘  =  -  0  0  e^*^^St)SCa,{t) 

t-#:(OQo(0  e-^^^Ci.t)SC^,i.t)  2A^ 

In  the  case  of  equal  chirp  rates,  a=  ^  ,we  have 

2SS^  ISSCe-'*^  -e'“'''"Qo(0 

2dSCe‘**  2dC"  0 

0  2A^ 

In  the  2-photon  resonance  case,  ^  =  0 ,  and  equal  chirp  rates,  a  =  ,  we  have 

0  0  -e'^'"Qo(0 

0  0  0 

0  2A^ 

We  can  see  that  a^(t)  =  a^(0) ,  and  we  are  left  with  the  system  of  2  differential 
equations  for  the  probability  amplitudes  aj^it)  and  h{t) .  Making  the  adiabatic 
elimination  of  the  excited  state  I T)  (assuming  that  hit)  «  0 ),  we  have 


and 


The  solution  is 


Therefore,  the  evolution  operator  for  the  2-photon  resonant  excitation  by  the 
equally  chirped  pulses  under  condition  of  the  adiabatic  elimination  of  the  trion 
state  has  the  following  form: 
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f 


U(0  =R 


0  0 

0  1  0 

0  0  1 

V  J 

e'*'  cos  i9  sin  i9  0 

e'*^  sin  &  cos  &  0 

0  0  1 

e'^cos"^>  +  sin"^i 
6’-'^^(e'^)-l)sin^icos^> 


R. 


0 

0 


j 


e  cos  i9 
-e'*^  sin  Q 
0 


0  0 
1  0 
0  1 

cos"^i  +  e'^  sin^^i 
0 


£,-'4 


e'^’’  cos  &  0 


g  -r  (g“ - 1)  sin  ,9  cos  i9  0 


0 


1 


(81) 


where  Sit)  =  f  Qo(r  )<ir  /  (4Ag) ,  and 

Jo 


For  the  qubit  states,  the  evolution  operator  can  be  written  as 


(  JS{t)l2 


U(0 


iS{t)l2 


COS  r9  +  ^ 


-iS{t)l2 


sin  r9 


^'^^(^'"(^)/'-^-*"(^)/")sinr9cosr9 


^-iS{t)l2 


cos^  3 sin^  3 


^  ^iSit)/2 


cos 


(5(0  /  2)  +  /sin(5(0  /  2)cos(2r9)  sin(5(0  /  2)sin (2*9) 


sin(5(0/2)sin(2*9) 


cos 


(5(0  /  2)  -  /sin (5(0  /  2)cos  (2*9) 


^^iS{t)/2  cos(5(0/2)  I  +  /sin(5(0/2)(cos(A^)sin(2*9)<T^ 


-sin 


(A^)sin(2*9) 


-Hcos 


_  ^iS{t)/2^-iS{t)  n-o/2 


(82) 


where 

n  =  (-cos  (  A^)  sin  (2.9) ,  sin  (  A^zi)  sin  (2.9) ,  -  cos  (2.9)). 

3.10  Electron  Spin  in  a  Quantum  Dot  as  a  Qubit 

In  the  previous  sections,  we  developed  several  methods  of  an  arbitrary 
manipulation  of  a  qubit  wave  function  using  the  geometric  phase.  Now  we  apply 
the  proposed  scheme  to  electron  spin  states  in  a  charged  quantum  dot.  Because  of 
quantum  confinement,  the  state  of  the  electron  can  be  expressed  as  a  product  of 
the  Bloch  function  and  an  envelope  function,  which  has  a  typical  scale  of  the 
quantum  dot  size,  a  few  nanometers.  The  energy  level  structure  and  optical 
selection  rules  have  been  discussed  in  the  literature.*^  A  commonly  accepted 
energy  level  structure  comprises  4  levels:  the  2  electron  spin  states  and  2  trion 
spin  states.  Figure  23  shows  2  arrangements  of  the  energy  levels  and  polarization 
selection  rules,  which  provide  a  possibility  of  optical  control  for  the  electron  spin 
qubit. 
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Fig.  23  Optical  selection  rules  in  different  bases:  (a)  the  mixed  basis  is  used,  where  the 
electron  spin  states  are  in  the  x  basis  while  the  trion  states  are  in  the  z  basis;  (b)  coupling 
scheme  is  in  the  x  basis. 

The  usual  experiment  of  the  electron  spin  control  is  performed  at  low  temperature 
(~  1  K).  An  external  magnetic  field  in  the  Voigt  configuration  (of  order  2h-7  T) 
is  applied  along  the  x  axis,  perpendicular  to  the  sample  growth  direction,  the  z 
axis.  Zeeman  splitting  of  the  electron  and  trion  spin  states  is  on  the  order  =  10 

meV  and  co^  =10//eV,  correspondingly.  At  these  conditions,  taking  into  account 

the  optical  selection  rules,  the  4-level  system  can  be  considered  as  a  double  A  - 
system.  This  coupling  scheme  is  shown  in  Fig.  23b,  where  we  indicated  by  H  and 
V  the  optical  field  couplings  with  the  orthogonal  polarization.  The  shown 
coupling  scheme  is  in  the  so-called  x  basis. 

An  alternative  arrangement  is  depicted  in  Fig.  23a.  In  this  case  the  mixed  basis  is 
used,  where  the  electron  spin  states  are  in  the  x  basis  while  the  trion  states  are  in 
the  z  basis.  Using  the  2  <T^  or  cr”  polarized  fields,  one  can  couple  the  electron 
spin  states  I  X_)  and  I ,  as  shown  in  the  Fig.  23a.  This  is  the  case  where  our 

3-level  model  can  be  implemented.  The  corresponding  Hamiltonian  has  the  form 
obtained  in  Eq.  48.  Assuming  large  detunings  of  the  pump  and  Stokes  field 
frequencies  from  the  transition  frequencies  to  the  trion  state  IT),  after  the 

adiabatic  elimination  of  the  trion  state,  for  the  case  of  completely  overlapped, 
identically  linearly  chirped  pump  and  Stokes  pulses,  we  obtain  the  following 
Hamiltonian 


h(  5 

-5  ^ 


(83) 
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where  S  =  (o^  -Aco ,  and  Aco  =  cOp-cOg,  =  0.pQ{t)Q.gQit){\  +  /  (2A) 

is  the  effective  2-photon  Rabi  frequency,  A  =  (Oj -cOp,  and  A<f>  =  ^p-^^.  Note  that 
here  we  have  used  //q^.  » i\j . 

Now  we  show  that  the  contribution  of  this  oscillating  term  can  be  neglected  when 
the  pulse  duration  is  longer  than  A<»“'  « .  To  demonstrate  this  and  justify  the 

procedure  of  the  adiabatic  elimination  of  the  trion  state,  we  numerically  solve  the 
time-dependent  Schrodinger  equation  with  the  Hamiltonian  in  Eq.  48.  That  is,  we 
compare  our  analytic  solution  with  the  exact  solution  of  the  Schrodinger  equation 
without  the  adiabatic  elimination  approximation.  An  example  of  the  comparison 
is  shown  in  Fig.  24. 


Fig.  24  The  population  dynamics  of  the  resonant  qubit  states  with  (dotted  lines)  and 
without  (solid  lines)  adiabatic  elimination  of  the  trion  state  in  the  3-level  system.  The 
excitation  is  generated  by  the  sequence  of  2  pairs  of  n  pulses  (Gaussian  pulse  envelops)  with 
the  relative  phase  =  nl2. 

The  example  in  Fig.  24  shows  the  dynamics  of  the  population  of  the  resonant 
qubit  excited  by  a  sequence  of  2  pairs  of  linearly  chirped  pulses  with  the  relative 
phase  between  the  pairs  being  ^  =  .tt  /  2 .  The  parameters  of  the  excitations  are  the 
transform-limited  pulse  duration  Tq  =  100  fs,  the  pulse  area  of  the  pump  and 
Stokes  pulses  is  equal  to  n,  liner  chirp  rate  a  /  Tq  =  20 ,  Atq  =  6 ,  =  1 ,  the 

ratio  between  the  maximum  Rabi  frequency  and  the  single-photon  detuning 
Qq  /  A  =  0.77 .  We  observe  a  negligibly  small  amount  of  population  in  the  trion 

state  I T)  at  the  intermediate  time.  The  presented  comparison  also  demonstrates  a 
reasonable  agreement  between  the  proposed  control  schemes  and  the  exact 
numerical  solution.  Note  that  the  total  time  for  the  qubit  operation  in  Fig.  24  is  on 
the  order  of  25  ps,  which  is  much  shorter  than  the  typical  lifetime  of  the  trion  state 
as  well  as  the  time  scale  of  other  forms  of  decoherence,  such  as  that  induced  by 
the  electron-phonon  interaction. 
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3.11  Possibility  to  Realize  Geometric  Gates  Using  Nitride 
Structures 

We  presented  the  analytic  expression  of  the  evolution  operator  of  the  electron  spin 
in  a  quantum  dot,  which  provides  a  clear  geometrical  interpretation  of  qubit 
dynamics.  Using  the  analytic  form  of  the  evolution  operator,  we  proposed  a  set  of 
single-qubit  rotations  that  are  solely  based  on  the  geometrical  phase.  Our  proposal 
combines  the  pulse  area  control  with  the  adiabaticity  by  using  chirped  pulses.  To 
estimate  the  time  scale  of  the  proposed  operation,  we  can  use  the  100-fs  pulses, 
which  become  picosecond  pulses  after  linearly  chirping.  That  amount  of  chirping 
is  sufficient  to  provide  adiabatic  excitation^^  and  can  be  readily  produced 
experimentally  using  commercially  available  laser  systems.  Using  parameters  of 
the  dipole  moments  //or,ri  ~  200  D  available  in  InGaN/GaN,^’  GaN/AlN^^’^^ 

quantum  dots,  and  detuning  A  =  5  meV,  we  estimate  the  peak  amplitude  of  the 
pulses  on  the  order  of  10^  -10^  V/m.  The  demonstrated  adiabatic  manipulation  of 
a  qubit  using  only  the  geometric  phase  has  some  advantages,  since  it  reduces  the 
requirements  of  perfect  tuning  of  the  control  field  parameters  and  is  significantly 
more  robust  against  noise. 

4.  Conclusions 


We  have  shown  that  we  can  explain  the  kinetics  of  the  photoexcited  carriers  in 
GaN  and  ZnO  by  including  only  carrier-carrier  and  carrier-phonon  interactions 
and  a  relatively  simple  2-band  electronic  structure  model.  Furthermore,  we  have 
presented  a  detailed  analysis  of  the  effect  of  different  dephasing  machanisms  on 
the  carrier  kinetics.  Finally,  we  have  experimentally  measured  and  simulated 
using  our  theoretical  model  the  temperature  and  pulse-power  dependence  of  the 
photoluminescence  spectra.  We  have  found  that  carrier-carrier  scattering  is  a 
principal  dephasing  mechanism  in  bulk  GaN  and  ZnO  for  excitations  below  the 
LO  phonon  energy  above  the  band  edge.  The  coherent  optical  control  of 
electronic  excitations  in  these  materials  will  require  confined  geometry,  e.g.,  a 
quantum  dot  heterostructure,  that  will  allow  us  to  reduce  the  available  density  of 
final  states  in  scattering  processes.  Our  theoretical  investigation  shows  the 
possibility  to  realize  geometric  quantum  gates  using  nitride  quantum  dots. 
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